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1. SFs RELEASES IN THE F REGION 

Sulfurhexafluoride (SF 6 ) is used primarily as a gaseous insulator in high-voltage systems. 

The insulation ability comes fhnn the fast removal of free electrons through electron attachment. 
SF 6 has been used as a release chemical in the earth's ionosphere to examine plasma instability 
processes. The chemistry of SFe in the ionosphere requires very few reactions for short time-scales 
due to the limited number of constituents at ionospheric altitudes. For longer time-scales the 
reaction list becomes complex. The region near the release becontes an electron hole and a 
negative-positive ion plasma. Earlier studies used simplified chemistry to model the electron hole 
that forms in short time [/{osenherg, \96A-, Mendilh and Forbes, \9%2\ Bernhardt, 1984]. 
Mendillo and Forbes [1982] assumed the major ion would be SFe', while Bernhardt [1984] 
predicted SFs' would dominate due to the high ambient temperature [Fehsenfeld, 1970], 
Experimental results from an SFe release in the E region suggested that F' may be the dominate 
negative ion shortly after the release [Hunton et al, 1987], 

The most recent ionospheric releases of SF« were done as part of the CURES (Combined 
Release and Radiation Effects Satellite) program. The NASA/Boston University rocket campaign 
was entitled CRRES-at-Kwajalein. Two rockets on separate evenings released SFs into the 
ionosphere to reduce conductivity through electron attachment and to mass-load a magnetic flux 
tube with heavy negative ions. The modified plasma produced a perturbation to the 
electrodynamics of the evening equatorial ionosphere where dramatic plasma irregularities develop 
naturally often nightly. The natural devel(q)ment of the irregularities in the equatorial F region has 
been termed Equatorial Spread F (ESF). Similar and separate releases of barium and HjO in the 
evening equatorial icmosphere were used in the Coloured Bubbles Experiment and the Brazilian 
Ionospheric Modification Experiment (BIME), respectively, to alter electrodynamics of the 
ionosphere [Haerendel et al, 1983]. Such studies of artificially-initiated ESF are expected to 
illuminate the natural causes of ESF. 

The purpose of this paper is to examine the current understanding of the SF 6 chemistry by 
comparing experimental data and modeling of two SFe releases from rockets launched from 
Kwajalein island during the CRRES-at-Kwajalein campaign. The two release rockets are named 
CRRES-1 and CRRES-2 within this document. The initial several minutes of the release are 
observed by airglow produced by the chemistry occurring within the release cloud. The late time 
diemistiy (IS minutes after the release) is observed by an Ion Mass Spectrometer (IMS) on a 
second instrument rocket. The combination of the two observations provides bracketing 
informatirm that will help discriminate between reaction sets and physical processes included in the 
model. 

The results of the study suggest that the collective chemistry presented by several authors 
\Mendillo and Forbes, 1982; Bernhardt, 1984, 1987; Hunton et al, 1987] correctly describes the 
^ort term chemical processes within the release. As time progresses, oxidation of the SF„ 
molecules increases the necessary reaction list. There are few laboratory studies on attachment 
characteristics of SOmF„ species but the electron density is not greatly affected by the presence of 
the large species list. For studies of plasma effects within the deep electron hole near the release 
point, SFe', SFs', and F' appears to be a sufficient list of negative ions. To examine late time 
chemistry and the broader spatial extent of the release one must include SOFs' and maybe SF 4 ' and 
SOF 5 -. 

This p^r presents the chemistry and airglow occurring within an SFe ionospheric release. In 
Section 2 the ambient background is modeled to provide correct context for modeling the release. 
Section 3 presents early time and late time observations and model results of the second SFg release 





in the CRRES-at-Kwajalein campaign. Section 4 provides a summary and advances suggestions 
for future work. 

2. SF# CHEMISTRY 

Mendillo and Forbes [1982] and Bernhardt [1984; 1987] review ion chemistry of ionospheric 
SFfi releases. Several important reactions are added by Hunton et al. [1987], It is apparent from 
previous and present SFe-ionospheric chemistry studies that there is still an element of unknown in 
the set of chemical reactions. The reactions used in this study are presented in Table 1. 

TABLE 1. SFfi Chemistry. 


DISSOCIATIVE ATTACHMENT 


SFg +e-^SFj- +F 

* 2 =V( 1+1 56x 10-3exp(4770/r„)] 



see Bernhardt [1984] 

SFg+e— ^F'+SFs 

^3 aO 

Mullen etai [1971] 

SOF4+e—+F 

00 

1 

0 

Sauers, communication 


RADIATIVE ATTACHMENT 

SFfi+e-^SFg- 

*a=(*,-^ 2 V 10 

see Bernhardt [1984] 

0 + e—^^0’ 

*5=1.3x10-15 

Branscomb [1964] 


CHARGE EXCHANGE 


0 ^ + SFg — ^ SFs"" + F + 0 

*6=1.5x10-0 

Fehsenfeld (1971) 

O++SF 5 —^5i->SF5^+0 

*7=10-0 

estimated (-k^) 

O2+SF5—^S->SF5^+02 

00 

0 

estimated (~k 5 ) 

SF 6 + 0 -^ 0 -+SF 6 

* 9 = 5 x 10-11 

Fe/»sert/eW[1971] 

SF5" + SFg *10 > SFe" + SF5 

*10=2.9x10-1'! 

Seppelt [1977] 

F'+SFe *11 >SF6~+F 

* 11 = 1 . 2 x 10-12 

Greenberg and Hargis [1989] 


CHARGE INTERCHANGE 


SFj-+ 0 —^^F +SOF4 

* 12 = 1 . 1 x 10-11 

Hunton et al. [1987] 

SFg" + SOF4 *13 > SOF5" + SF5 

*13=8.5x10-0 

Van Brunt el al. [1988] 

SFg" + SF4 SFj- + SF5 

*14=3.7x10-10 

Babcock and Streit [1981] 


ION-ION RECOMBINATION 

0++0- -*I 5 >o*+o 

*15=2.7x10-2(T/300)-5 

Olson et al. [1971] 

O^+SFe—^aii-> 0 *+SF 6 

*16=6x10-*(T/300)-5 

Bortneretai. [1972] 

SFs"’ + SFg -*-l7-> SF5 + SFg 

*17=3.9x10-8(T/300)-5 

Church and Smith [1977] 

X++Y~ *<* >X + Y 

*18=6x10-*(T/300)-5 

Bortneretai. [1972] 


E-ION RECOMBINATION 


SF 5 '*^+j -^ -■»SF 4 +F 

*ig=2xlO-2(T/30O)-5 

Biondi [1973] 


ASSOCIATIONAL DETACHMENT 

0-+0_*2Q.^02+e 

* 20 = 2 x 10-10 

Fehsenfeld et al. [1969] 
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1 4BLE 1. SF6 Chemistry, (continued) 


COLLISIONAL DETACHMENT 


SFg'+M- 

-^SU->SF6 + e+M 

*21=1x10-12 

Wangetal. [1989] 

SFs'+M- 

-^2i-^SF5+e + M 

*22=1x10-*^ 

Wemgetal. [1989] 

F'+M— 

^^22_>F + e + M 

*23=5x10-1* 

Wangetal. [1989] 



COLLISIONAL DISSOCIATION 

SFfi'+M- 

-^a^SFj+F + M 

*24=1x10-13 

Wang et al. [1989] 

SFfi'+M- 

-^aS-^F’+SFs+M 

*25=1x10-13 

Wangetal. [1989] 

SF 5 +M- 

-^a^->F-+SF 4 +M 

*26=2x10-1'* 

Wangetal. [1989] 



PHOTODETACHMENT 


SF/-iai 

— >SF4 +e 

*27=1.0 

Bernhardt [1987] 

SFj-—^ 

—> SF 5 + e 

*28=0.5 

Bernhardt [1987] 

SFfi-— ^ 

->SF 6 +e 

*29=2.0 

Bernhardt [1987] 

0-_J2fi_ 

>0 + e 

*30=1.15 

Burch etal. [1959] 



NEUTRAL REARRANGEMENT 

O + SF 5 - 

^ai-^SOF4+F 

*31=2x10-11 

Van Brunt [1990] 

SF 5 +SF 5 

-^22->SF4+SF6 

*32=5x10-13 

Van Brunt [1990] 


2.1 Electron Attachment 

Initially, the dominant reaction is SF6-e attachment, which proceeds at a rate kj = 2.2xl0~^ 
cm^s*^ for low energy electrons (Tg-0.1 eV) [Fehsenfeld, 1970; Mothes et al., 1972; Spence and 
Schultz, \9Ti, Phelps and Van Brunt, 1988] 

SFg + e—Products 

The cross section dramatically drops off with increasing electron energies [Mullen et al., 1971; 
Kline et al., 1979; Phelps and Van Brunt, 1988]; only the lowest energy electrons have a 
significant attachment cross section with SFe. In a release in the nighttime ionosphere, SFe 
vibrational modes can significantly cool electrons to provide low energy electrons for rapid 
attachment and k| can be considered constant. The products of the electron attachment reaction 
depend upon the vibrational state of the parent SFe molecule prior to attachment. The branching of 
SFe attachment chemistry can be modeled with a ground state SFe' and a vibrationally excited SFe* 
*. SFe'* represents all SFe* molecules with vibrational energy below the F-SFs* dissociation 
energy. Bernhardt [1984] uses three vibratiorudly excited species, SFe* , SFe* and SFe* , but 
the two additirmal excited species dissociate almost instantaneously into end products. 

The three main branches of the attachment process are 


SFj+e- 

-^SF,* 

(2.1) 

SFg +e- 

-^>->sf;+f 

(2.2) 

SFj +e- 

-^F+SF, 

(2.3) 
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where ki = k^. The branching of the attachment process is temperature dependent. 

Reaction (2.3) is negligible at ionospheric temperatures. The temperature dependent rates of 
reactions ( 2 . 1 - 2 ) are i^proximated with 


k^ = +136 X 10-^ exp(4770 / T„)] 

k^ = k, - k 2 = ^,/[l + 640exp(-4770/T„)] 


(2.4) 


(2.5) 


[Bernhardt, 1984; Fehsenfeld, 1970]. The metastable molecule SFe'* can auto-detach, radiatively 
stabilize, or collisionally stabilize, as illustrated in the following reactions. 


SFj*-^!^SF,+e 


( 2 . 6 ) 


SFg*—^SFj+hv 


(2.7) 


sf;*+m—^^SF;+M 


( 2 . 8 ) 


Auto-detachment of SF 6 ~* is rapid: = l/a^ = lOps to 93ps depending on SFe vibrational energy 
[Odom et al. , 1975; Foster and Beauchamp, 1975; Lifshitz, 1983], CoUisional stabilization has 
a rate of k^ = 5.4x10''° ernes'*. Above 125 km altitude collisional stabilization can be n^lected. 
Radiative stabilization is approximately ten times slower than auto-detachment [Foster and 
Beauchamp, 1975], Thus, the net production of SFe* through the combined reactions of (2.1,2.6- 
2 . 8 ) becomes 


SF« +e—^^SF; 


‘6 
where 


(gr+MM} . 


(2.9) 


( 2 . 10 ) 


The electron attachment reactions dominate the very early time chemistry that defines the negative¬ 
positive ion plasma within the cloud. 

In the CRRES-at-Kwajalein SFe release, the IMS possibly observed SF 4 ~ at about the same 
density levels as SFf' and SFs~. Previous studies of SFe releases did not allow for SF4~ production 
>^ch can arise primarily through electron attachment to SF4. SF4 attachment chemistry is similar 
to SFe chemistry with attachment branches 


SF 4 +e-►SF 4 -' 


( 2 . 11 ) 


SF 4 +e-►F' +SF 3 (2.12) 

The attachment rate for SF4 11 times slower than SFe attachment (Af/7/er et al., 1994). As with 
SFe radiative stabilization must occur. Auto-detachment rate is approximately 100 times fiister 
than radiative stabilization which make SF4* production very slow in a ratified atmosphere (Sauer 
etal, 1993). 

The apparent IMS observation of SF4~ (108 amu) densities may instead be SOFs' (105 amu) 
due to the relatively low mass resolution of the instrument. Though SFe* and SFs* will be 
important during the first several minutes of the release, other neutral molecules appear through the 
oxidation of SFn and more exotic negative ions emerge. In particular, modeling suggests that 






S 0 F 3 ' is more dense tiian SFe* and SFj* 10 minutes after the release. The main reaction producing 
SOFs" is 


SOFj+e->SOFj+F (2.13) 

2.2 Charge Transfer and Charge Interchange 

The instrumented rocket was launched approximately IS minutes after the release rocket 
launch. IS minutes after the SF 6 release reactions slower than the attachment processes have 
altered the ion constituents significantly. In particular, charge exchange will have altered the 
positive component within the release cloud through 

+ SFg-►SFj* + 0 + F (2.14) 

O* + SFj-►SF/ + 0 (2. IS) 

These are potentially important reactions for a daytime release where photo-detachment reduces the 
negative ion densities. The introduction of positive molecular ions will alter the electrody nam ics 
through mass loading of the field lines as well as plasma reductions through dissociative 
recOTibination. 

Two major adjustments to the native ion populations result from the dominance of 0 in the 
upper atmosphere. 

SFj+0->F'+SOF^ (2.16) 

SFj'+O-►O'+SFj (2.17) 

Reaction (2.3) is particularly important for the late tune observations of the negative ion plasma. 
The dominant negative ion at late times is F* with reaction (2.3) as the primary source of F'. 

Other n^ative ion species appear through fast interchange reactions as oxidation of SFq 
continues: 

SFj + SOF 4 -> SOFj + SFj (2.18) 

[Van Brunt et ai, 1988]. This leads to a potentially complex n^ative-positive ion plasma. 

2.3 Recombination 

As the plasma is altered from an atomic oxygoi ion-electron plasma to a molecular ion 
plasma, the recombination rates are accelerated. The O'^-e recombination is so slow as to be 
negligible (wlO*^^ cm*^ s'^). As positive molecular ions appear the plasma recombination rate 
increases to about lO'*^ cm"^ s*^ [Biondi, 1973]. Ion-ion recombination is not much slower at 
about 10** cm*^ s*^ [Bortner et ai, 1972]. 

2.4 CID and Other Reactions 

As time progresses additional loss reactions can enter into the chemistry. Hunton et ai 
[1987] suggest coUisionally induced detachment and dissociation as loss mechanisms for negative 
ions. Recent investigations [ITong^ro/., \9%%Olthoffet ai, 1989] provide the cross sectional 
information to allow estimates for rates at atmospheric temperatures (T~1000K). However, these 
reactimis are slow (< 10*12 cm^s'l), therefore, they probably have limited influence on the 
chemical perturbation except during late times after the release. 

Associational attachment generally are much fester. For example. 
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(2.19) 


O' +0-+e 

at 2x10" 19 cm'^ s'l [Fehsenfeld et al., 1966; 1969] and 


F'+H->HF + e (2.20) 

at 1.6x10*^ cm*^ s*l [Fehsenfeld et al., 1975] may be (grating with other negative ions though no 
rates could be found in the literature. Is it possible that the following are important? 


F'-t-O- 

-^OF+e 

(2.21) 

F -i-F- 

--—►F 2 +e 

(2.22) 

SF, +0- 

-^SOFj +e 

(2.23) 

SFj+F- 

-L^SF«+e 

(2.24) 


Because the rates are possibly near gas kinetic, thQ' could influence negative ion densities. 

2.5 Neutral Rearrangement 

Potentially important neutral reactions are missing in the review papers of ionospheric release 
of SF6. The neutral reactions may influence the late-time development of an ionospheric 
peituibaticn. Due to the dominance of 0 at F re^on altitudes, the oxidation of SFp molecules 
leads to several SOmFo-i neutrals [Van Brunt, 1985; Van Brunt andJ.T. Herron, 1990]. In 
addition, the ion-ion recombination, as well as the charge interchange chemistry given above, 
produces marry contaminant neutral species withui the chemically perturbed r^on. Some 
contaminarrt species are electron attachers (SF4, SFs, SO2F, SO2F2, SOF4). 

2.6 Two Regions of Chemistry 

Early after the release, the SF« is cool (T~300K) and dense (nspe > n,). Before the SFe release 
cloud expands to a 20 km radius, the rapid attachment chemistry has occurred and a deep electron 
hole has been created. The resulting ion-ion plasma cloud may increase the mass or decrease the 
(xmductivity on the magnetic flux tube to initiate instability processes at the edge of the plasma 
doud. Additionally, the recombination of SF«~ and O'*' has oiough energy to create highly excited 
O and lead to visible aitglow emissions. The important chemical processes are illustrated in Figure 
1 . 

In the broader r^on after the release cloud has expanded and collisionally heated, the 
diemistry leads to a diffuse negative ion perturbation that does not affect the electrodynamics. The 
recombination of SFs* or F* with 0*^ does not lead to visible aitglow emissions but does produce 
UV emissi<Mi. The important chemical processes are illustrated in Figure 2. 
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0*+ hv(7774A, 8446A) 
O +hv(1304A, 1356A) 


If SFg is cool and dense, then: I. SF^, 0‘, O 

n. Deep electron hole. 

III. Visible and UV airglow. 


Figure 1. Chemistry for a dense cool SFe release in the Fj region. 


XX. 

SF/ sf; 

XX X 

F' O 

4 . ^ 

X X. 

X 


SOF^.F- 


neutrals 


0 

O* 


O +hv(1304A, 1356A) 


+ r^ + 


If SFg is hot and less dense, then: I. SFj,F’, SF^O 


II. Broader region disturbance. 

III. UV airglow. 


Figure 2. Chemistry after a SF« release has expanded and heated. 
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2.7. SFe-HjO Competition 

The release rocket and observation rocket outgas atmospheric gases that are imbedded in the 
skin of the rockets. The primary constituent of the outgassing is H 2 O. The HjO competes with 
SFe in die positive-ion chemistry. HjO and charge exchange and subsequent hydration 


chemistiy occur. 

HjO+0^-►HjO^+O (2.25) 

H jO+H jO^-►HjO* + OH (2 26) 

The increase of molecular ions leads to increased recombination of the plasma. 

HjO*+e-►OH + H (2.27) 

HjO^+e->HjO + H (2.28) 

Within the SFe induced perturbation the electron density is reduced and ion-ion mutual 
neutralization will be the avenue for H 2 O''' and HjO'*' loss. 

HjO*+M--^HjO + M (2.29) 

H,0^+M--►HjO + H + M (2.30) 

A further complication of the SF 6 release chemistry is the possibility of water clustering of the 
negative ions. 

SF„-+HjO + M->SF„-(HjO) + M (2.31) 


However, a third-body is required for hydration to proceed. Hydration will only be important for 
releases below 100 km and in the first moments after a release. 

The introduction of H 2 O will affect SFe chemistry through competition for ions and 
through indirect reaction of H 2 O chemistry products (OH, H) with SF„ fragments. Effects of H 2 O 
in an F region release of SFe is small. The only observable change is the increased presence of 
H 20 '^ as a minor i<Hi species. 

2.8. Airglow Emissions from an SFe Release 

Airglow observations were made during the CRRES-at-Kwajalein campaign. An all sky 
camera of the Center for Space Physics of Boston University was used to observe the spatial 
distributiem of the expected airglow enhancement during the SFe release. The enhancement is 


primarily fhrni the ion-ion recombination 

SFj- + 0*-►SFs + 0 + 13.06eV 

(2.32) 

sf;+o*- 

—►SF 4+0 + 12.1 lev 

(2.33) 

F+O* — 

-+F+0 + 10.20eV 

(2.34) 

SFj-+0*- 

-^SF 5 + 0 + 9.91eV 

(2.35) 

S 0 F 5 ^+ 0 * 

-►SOFj+O + S^leV 

(2.36) 
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Recombination of with SFr and SFe* can produce highly excited states of atomic oxygen. 
These excited states decay in multiple steps. The specific excited states of interest and their 
associated emission lines are; 


0(3p*P)- 

—>0(3s*S) + hv(7774>!\) 

(2.37) 

0(4p*P)- 

—>0(3s’S) + hv(4368A) 

(2.38) 

0(3p*P)- > 0(3s^S) + h v(8446;\) 

The lower states relax to ground state via 

(2.39) 

0(3s*S) — 

->0(3p^P) + hv(1356A 

(2.40) 

0(3s’S)— 

-^O(2p^P) + hv(1304;^) 

(2.41) 


The branching ratios for the excited states of O produced in the recombination have been 
approximated by using the O'^'-e recombination branching ratios [Bernhardt et al ., 1986]. 

Recombination of O with SFs* and F~ can only produce the lower excited states of atomic 
oxygen diat radiate ultraviolet airglow. Recombination with SOF 3 ' can only produce the low 
metastable states of 0(^5) and 0(^0). The differing energetics of the recombination reactions 
allows for interesting spatial structure of the airglow in the different wavelengths. In the small 
region w^iere SFs' densities are high, 7774A is bright, but in the broader region where SF5* and F~ 
dcHninale, 13S6A is important. It would be instructive to measure 7774A and 8446A or 1304A 
and I356A lines together. 

3. MODELING 

The NASA/Boston University CRRES-at-Kwajalein campaign took place during the summer 
of 1990. Two sounding rockets (CRRES-1 and CRkES-2) released SFe on separate evenings in 
an attempt to induce equatorial spread F phenomena. Multiple diagnostics were used to examine 
the resulting perturbations to the ionosphere including Boston University’s all sky optical camera 
and incdierent scatter radar (ALTAIR). An ion>mass spectrometer (IMS) from Phillips 
Laboratory Directorate of Geophysics flew on an instrumented rocket launched fifteen minutes 
after the second release rocket (CRRES-2). Table 2 compiles the geographic, solar, and 
atnrospheric ccMulitions for the equatorial launches. Both flights released approximately 35 kg of 
SFe (1.4x 10 “ molecules) at a temperature of 320K. The release canister took 3 to 7 s to empty. 

The resulting SFe-ionospheric chemistry creates four categories of alterations to the ambient 
conditions: 

i. Electron depletion. This interesting perturbation alters electromagnetic wave 
scattering and propagation. 

ii. Neaative-positive ion plasma. The replacement of electrons by heavy negative 
i(xis alters the transport properties along magnetic field-lines, increases cross 
field conductivity, and increases the gravitationally-induced current. 

iii. Plasma depletion through ion-ion recombination. The plasma reduction takes 
place on the time scale of minutes and eventually decreases conductivities. 

iv. Neutral composition perturbation. There are several metastable constituents 
that participate in airglow enhancements. 









The study of the SFe release cannot be done exclusively with models of the chemical 
perturbation. Comparismi of simulation and experiment would have only limited validity. To 
obtain more quantitative comparisons, the dynamics and electrodynamics of the CRR£S-at- 
Kwajalein release were modeled iit conjunction with the chemistry. 

TABLE 2. CRRES*at-Kwajalein release conditions. 



CRRES-1 

CRRES-2 

Launch date 

11 August 1990 

14 August 1990 

Launch time 

08; 10:00 UT 

08:29.00 UT 

Release time 

08:17:20 UT 

08:34:34 UT 

Release Altitude 

375.5 km 

428 km 

10.7 cm Flux 

187 

200 

<10.7 cm Flux> 

228 

228 

Sunspot Number (Rj|) 

161 

229 

Geo-Longitude 

166.90“ E 

167.0898“ E 

Geo-Latitude 

10.97“ N 

10.8825“ N 

Magnetic Dip 

~0“ 

~0“ 

Ap Index 

9 

26 


3.1 Ambient Conditions 

The atmospheric conditions used as a background for the modeling studies were determined 
by MSIS90 [Hedin, 1991]. The pre-flight electron densities and seasonally averaged plasma 
velocities were obtained from Altair radar measurements (Peter Sultan, personal communication) 
(Figure 3). 



Electron Density (cm-3) Local Time (hours) 

Figure 3. (a) Electron density profile for 15 August 1990 measured between 8:23;24-8;28;40 UT, 
about 9 minutes before launch of the release rocket, (b) Plasma drift velocity averages determined 
from Altair data taken during summer 1990. 
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The IMS measurement indicates negative ions concentrations are small compared to the total 
positive ion concentration. Thus, the IMS^neasured O'**, NO'*', O 2 *' profiles should app.'-oximately 
represoit the undisturbed ccmditions. Unfortunately, the IMS chaimel was saturated and does 
not provide information on the '^e" profile above 300 km. The IMS data combir^ with electron 
profoes obtained fron AHair Radar data do give a consistmt picture of the ionosphere at the time 
of the instrument rocket fly-4hrough. The electron profiles for 19.57 LT (8.43 UT) and 19.97 LT 
(8.83 UT) and the IMS measurements at s^rproximately 19.97 LT (8.83 UT) suggest that the 
ionosphere dropped at an average rate of IS m/s (Figure 4). The seasonal average vertical velocity 
is upward during the experiment time interval. Perlu^>s the seasonally averaged vertical velocity 
does not accurately represent the vertical velocity for the night of the release. This is not surprising 
due to the extreme variability of the vertical velocity near sunset. 

We attempted to model the ambient ionosphere fimn the time of the Altair measurement to the 
time of the IMS measuretiKnt. Comparison of model results to positive ion chemistry 
measurements will provide additicmal insight into the correctness of the ionospheric vertical 
velocity history. Three simulations were performed to reinforce the thesis th^ the ionosphere was 
moving downward ~15 m/s after the SF6 release of CRRES-2. The simulations use; the seasonal 
averaged vertical velocity (Figure 3), zero vertical velocity and, constant downward velocity of • 

IS m/s. The simulations b^in with conditions provided by MSIS90 neutral atmosphere and the 
electron profile measured by the Altair radar at 19.57 LT. The simulations continue until the time 
of the ot:^rvation rocket apogee (19.95 LT or 900 s after the SF6 release). The first simulation 
that uses the seasonal average vertical velocity produces a very high F regirm ledge and a high 
NO*^ profile "knee" (Figure 5a). Using no vertical motitm, the F rqpon ledge and the NO'*' knee 
are better but still the electron and O'*' profiles are too high (Figure 5b). Using vertical velocity of 
-IS m/s the simulation and experiment data have a good match in total electron density (Figure 5c). 
The Oa*** and NO'*' density profiles are approximately correct in sh^ but the lower altitude 
magnitudes do not compare fiivorably wifti the data. The total molecular ion density observed is 
approximately 3 times the modeled doisity. This could be due to differences in the true conditions 

frcHn the model conditions (for example, T«, [N 2 ], [O ]). 



Figure 4. The IMS ion data and the Altair radar electron profiles suggest that the ionosphere 
dropped during the CRRES-2 release experiment at an average velocity of -IS m/s. The O'*' ion 
data saturated but the F region ledge in foe later electron profile and foe O'** data are approximately 
20 km below foe earlier ledge. 
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Seasonal Averaged Vz 




Figure S. Simulations of the ambient chemistiy of the CRRES-2 release experiment. The 
simulations all began with same initial conditions. The initial electron profile is identical to the 
Altair measurement in Figure (3). (a) The seasonal average vertical velocity is used. TheF 
region ledge and the NO'*' "knee” are too high, (b) With no vertical motion the F region ledge is 
still too high, but the NO*^ "knee" appears at the correct altitude, (c) With a constant vertical 
velocity of *15 m/s the simulation and experiment data have the best match. 


12 












With a reascaiable estimate of the electron profile and vertical velocity time history one can 
b^in to compare the model results of the SF« release with observational data. The constant 
downward plasma velocity of-IS m/s will be used in the following studies. This velocity 
reproduced die observed ambient electron profile at the time of the observation rocket fly-through. 

3.2 Motion of the CRRES SFt Release 

To properly compare model results with airglow and IMS density measurements one must 
determine tte path of the neutral and ion clouds created by the release experiment. This requires 
knowledge of positi(m and velocity of both the release rocket and the instrument rocket. Table 3 
lists trajectory information for the release rocket against mission elapsed time (MET). Table 4 lists 
the instrument rocket information. 


TABLE 3. Release rocket trajectory against mission elapsed time (MET). 


MET(s) 

Latitude 

Longitude 

Altitude 

0.00 

9.4012 

167.4651 

0.037 

2.40 

9.3998 

167.4570 

2.017 

30.00 

9.4652 

167.4441 

37.019 

60.00 

9.6148 

167.4150 

112.676 

90.00 

9.7611 

167.3841 

179.848 

120.00 

9.9044 

167.3511 

238.722 

ISO.OO 

10.0459 

167.3169 

289.423 

180.00 

10.1853 

167.2815 

332.104 

210.00 

10.3228 

167.2455 

366.883 

240.00 

10.4596 

167.2084 

393.786 

270.00 

10.5953 

167.1705 

412.914 

300.00 

10.7301 

167.1331 

424.334 

330.00 

10.8669 

167.0963 

427.334 

334.00 

10.8825 

167.0898 

427.942 

360.00 

11.0038 

167.0596 

424.191 

390.00 

11.1406 

167.0228 

412.582 

420.00 

11.2777 

166.9861 

393.237 

450.00 

11.4143 

166.9493 

366.310 

480.00 

11.5511 

166.9125 

332.639 

510.00 

11.6880 

166.8758 

289.040 

540.00 

11.8248 

166.8390 

238.399 

570.00 

11.9617 

166.8023 

179.616 

600.00 

12.0985 

166.7655 

112.502 

630.00 

12.2353 

166.7287 

37.284 

660.00 

12.3722 

166.6920 

2.767 
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TABLE 4. Instniment rocket infomution against mission elapsed time (MET). 


MET(s) 

Latitude 

Longitude 

Altitude 

0.00 

9.4012 

167.4651 

0.037 

S.OO 

9.4029 

167.4663 

2.017 

30.00 

9.4366 

167.4974 

36.552 

60.00 

9.5163 

167.5706 

111.488 

90.00 

9.5944 

167.6396 

177.946 

120.00 

9.6709 

167.7050 

236.134 

ISO.OO 

9.7460 

167.7675 

286.192 

180.00 

9.8201 

167.8275 

328.234 

210.00 

9.8933 

167.8855 

362.365 

240.00 

9.9656 

167.9419 

388.668 

270.00 

10.0375 

167.9972 

407.196 

300.00 

10.1089 

168.0516 

417.999 

327.00 

10.1732 

168.1002 

421.130 

360.00 

10.2515 

168.1601 

416.491 

390.00 

10.3228 

168.2150 

404.182 

420.00 

10.3944 

168.2709 

384.137 

450.00 

10.4664 

168.3280 

356.310 

480.00 

10.5391 

168.3870 

320.639 

510.00 

10.6125 

168.4483 

277.040 

540.00 

10.6871 

168.5124 

225.399 

570.00 

10.7628 

168.5798 

165.616 

600.00 

10.8400 

168.6512 

97.502 

630.00 

10.9097 

168.7286 

26.284 

660.00 

10.9063 

168.7390 

15.767 


A simple and computationally &st release cloud model was developed to predict approximate 
position, d^ity and size of a neutral release cloud (ignoring chemistry). The 427.95 Im altitude 
is above the r^txi \^re difiusive expansion through the background gas can be used to 
accurately describe the early expansion and titll. Therefore, a model wdiich included momentum 
terms in the fidl and expansion was necessary to model accurately. It is assumed that diffusion is 
inportant and a Gaussian cloud is used rather than a shell cloud. A center-of-mass (COM) model 
was developed vidiich follows the motion of the center of the (jaussian release cloud as it falls 
through an MSIS defined atmosphere. The expansion and heating of the cloud are also computed. 
The basic assumption is that the cloud maintains a Gaussian shape described by 


n(x,y,z,t) = M„(Ocxp 




(3.1) 


wdiere x, y. z represent eastward, northward and vertical directions, n is density, and is peak 
doisity. The p^ density goes as 


«a(0 = 


N 


r +r » 


(3.2) 




■x-y 2 

The cloud center position (Xo,yo.z^ is followed with simple center-of-mass equations of motion 
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(3.3) 



Uo is Ae cloud center velocity, K, is the collision frequency between the release gas and the 
badcground neutrals, and g is gravitational acceleration. The expansion of the cloud is determined 
for each limb of die cloud. The horizontal limbs (r^ and ry) are assumed to expand identically. 

The upper and lower limbs of the cloud, and r^., are assumed to have independent expansion 
rates. The velocity of eiqiansion is determined by the velocity of the n^e radius of each cloud limb. 
This radius grows at the rate defined by the velocity of the gas at that radius, which can be 
determined by the momoitum equation 







(3.4) 




= -g- ('■,+ ) - (^O + 'i. )] - ■ 


* 


+v)l- 


where Uf^ is velocity of the horizontal limb expansion, Uf 2 + is the velocity of the upper limb 
expansion, is the velocity of the lower limb expansion, and m, is the release gas mass, and To 

is the release gas temperature. The pressure term in each equation reduces as the Gaussian limb 
radii in the denominator increase. The collisional heating of the cloud can also be approximated 
throu^ 


c, dt 






3^(7;-7;)-^(u„-uJ^ 




(3.7) 


Only the center of the cloud is used to define the tonperature of the vdiole release cloud. The 
solutions assume that the background gases are unchanged by the release. This is true for the large 
scale problem. Only at early times (t<10s) will the release affect the background densities and 
temperatures. The CRRES-2 SF« release is modeled using the above COM model with the 
MSIS90 neutral atmosphere model. The northward and eastward rocket velocities at release were 
503 m/s and -136 m/s, respectively. 

Figure 6 gives the expansion history and temperature of the modeled release cloud. The 
temperature increase is slow due to the internal energy capacitance of SFe. The slow increase 
allows SF6' to be produced as the dominant ion for a significant period. The cloud moves initially 
with the rocket trajectory and expands adiabatically for one mean free path (~ SO km). After two 
minutes, dififosion controls the motion. Within five minutes, the release cloud has the same 
horizontal motion and temperature as the background atmosphere. The region of chemical 
influence of the cloud has a radius greater than 100 km. The differential neutral drag force on the 
tq> and bottom limbs of the cloud causes the cloud to flatten. 
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Figure 6. (a) The expansion of the SF6 cloud limbs, (b) The coilisional heating of the SFe release 
cloud, (c) The peak density of the release cloud. 


The dense region of the ion cloud should closely follow the 7774A emission region measured 
by the Boston University all-sky c<unera. The horizontal posititm of the emission region center 
versus time after release of Table 5 was obtained from the a!! sky camera images (Michael 
Mendillo, personal communication). The altitude column of Table S is the predicted altitude of the 
ion cloud center given an downward plasma drift of 15 m/s. The North-South motion of the 
7774 A emission should approximate the neutral gas motion since this is parallel to the magnetic 
field. The East-West motion of the plasma perturbation should couple very rapidly to the 
bacl^round plasma motion. It is assumed that the ion cloud is formed within a chemical time 
omstant of ^ attachmoit process (10s). The background plasma East-West motion should be 
approximately 80% of the neutral wind. Additionally, the vertical motion of the ion cloud is 
assumed to be coupled to the background plasma motion. Altair measurements suggest that there 
was no vertical motion of the ionosphere for up to an hour after the rel^e. The background 
plasma was moving approximately 130 m/s East and 135 m/s South. These numbers were 
obtained fixxn the observed motion of the 7774A emission region. 


TABLE 5. Airglow center (Michael Mendillo, personal communication). 


TAR(s) 

Latitude 

Longitude 

Altitude (km) 

42 

11.18 

167.21 

426.1 

68 

11.30 

167.24 

425.7 

94 

11.34 

167.28 

425.4 

120 

11.38 

167.28 

424.4 

146 

11.34 

167.31 

424.6 

172 

11.26 

167.35 

424.2 

198 

11.22 

167.38 

423.8 

224 

11.22 

167.38 

423.4 

250 

11.22 

167.45 

423.0 

276 

11.14 

167.45 

422.6 

302 

11.18 

167.49 

422.2 


16 












The relative positions of the rockets and perturbation clouds are shown in Figure 7. The 
predicted ion cloud positi<xis produced by the COM model have an ofl^ of about 17 km east and 
17 km Dortii from the 7774A measurement positirms. The Boston University all sky camera 
position contains uncertainty of approximately 10 km (M. Mendillo, perstmal communication), 
which does not account for the entire ofi^. If tiie camera position is correct, the velocity required 
for the gas to reach the first 7774A measurement position is about 600 m/s, which is larger than 
the rdease rocket velocity. The only way this apparent problem may be overcome is if the gas 
release was directed, i.e., if the release takes place from a single no^e into space. The gas would 
accelerate out of the release tank to obtain a Erected velocity of 180 m/s (thermal velocity at 
300K) wifo respect to the release rocket velocity. If the directed velocity out of the release tank is 
northward tiien the COM prediction could better approximate what was observed. 

Another COM model run was done using a directed release of 130 m/s North and 130 m/s 
East witii respect to the release rocket. Figure 8 shows the foilure of this particular directed 
release to account for the 7774A measured positions. However, the later time positions obtained 
from the model are closer to the measured positions. 

Using the instrument rocket path and the model predictions of the ion cloud path, we note that 
the rocket-to-icm cloud separation varies only slightly betweoi these cases. The path of the ion 
cloud differs primarily in the initial minutes after the release. Thus, only slight variations are seen 
in the ion profiles predicted by several possible ion cloud paths. This provides confidence in the 
comparisons of data and model results. 

The e}q)ansion after the release may allow for freezing of the SF6 gas. The thermal energy of 
the gas is transferred to directed acceleration from the release rocket. This produces clustering or 
freeziitg of the SFe molecules and reduces the 7774A emission signature from the experiment. 


CRRES-2 



Longitude 


Figure 7. The position of release rocket, instrument rocket, ion cloud and neutral cloud as 
predicted by the COM model. The circle is the release point. The triangles are the measured 
positions of the 7774A emission in Table S. 
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Figure 8. A directed release is modeled with CCM program. The ion cloud position is still not 
able to match the 7774A emission positions from Table 5. 


3.3 Chemical Perturbation in CRRES-2 

Due to the quick coupling of the cloud to the neutral background and the lack of any 
significant perturbed plasma motitm, a simplified three-dimensional model of chemistry and 
dynamics can be used to simulate the diemical perturbation of the SFc release. Using axially 
symmetric cylindrical geometry in the frame of the ambient neutral wind, the three-dimensional 
aspects of the neutral motion can be reduced to two-dimensirms. Because the plasma is bound to 
the magnetic field lines, its only motion is parallel to the magnetic field. This caimot be properly 
model^ in the cylindrical geometry. However, the plasma perturbation can be modeled 
approximately by solving plasma motion in the plane containing the magnetic field and altitude 
directions witii rectangular coordinates. The combined cylindrical-rectangular, two-dirrwnsional 
model effectively rqrroduces a three-dimensional transport result. This simplified transport allows 
one to use large di^stry sets with less computational demand. 

The set of reactions used in the two-dimensional study is larger than in previous studies of 
irmospheric releases of SF«. The reactions of Table 1 are included with several ambioit ion 
reactions to develqr airglow predictions and density profiles. Thermodynamic calculations 
indic ate<< that 40% of the SF 6 may be frozen upon release. The calculation involved large 
uncertainties. For this reason, three cases were run using 100%, 10%, and 1% fiee-vapor of SF 6 in 
release. 

The tTKxlel results were compared to airglow measurements in East-West cross sections, and 
North-South cross sections of the 7774A airglow. The modeled 7774A airglow for North-South 
and East-West cross sectimis are plotted for the times of several airglow data points — 16,68, 120, 
and 250 secrmds after the release. Two model rruis of the release are presented in Figure 9 and 10 . 
Figure 9 diows the very low emission rates of tlra SF« release chemistry if the SFe is initialized 
wiA tire ambient upper atmosphere temperature of about lOOOK. This was suggested by 
Bernhardt et al. [1984]. Hot SFe leads almost entirely to SF 5 ' ions. The excited oxygen atoms 
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produced from SFs~ and O'*' iecombinati<Hi produce ultraviolet emissicMis. The atmosphere is 
optically thidc to diese UV emissions and could not be observed from the ground. Figure 10 shows 
tte very high emission rates of the SF< release chemistry. The SFe is initialized at the release 
temperature of300K and allowed to heat through collisions with the background gases. This 
allows SF6* to be produced for a significant time period (see Figure 6b). The 7774A intensities in 
a cold release of SFe are approximately 7 and 1.3 times larger than observed at 16 s and 250 s, 
respectively. 
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Figure 9. East«West and North>South cross sections of the 7774A airglow from the SF^ release 
cloud. The 7774A intensity is only 4 Rayleighs above the background for this model run. The 
main characteristic of the run is the assumed initial high temperature of the SFe gas. The SFe 
tenqwiature is set at the ambient atmospheric temperature (~1000K). 
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Figure 10. East-West and North-South cross sections of the 7774A airglow from the SFe release 
cloud. The 7774A intensity is 800 Rayleighs above the background in this model run. The main 
characteristic of this model is the assumed initial low temperature of the SFe gas. The SFe 
temperature is set at the release temperature (~300K) and allowed to collisionally heat in the 
model. 
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It may be su£Scieot to describe present observations if some thermodynamic considerations 
are taken into account. An ^>parent problem with interpretation of Red Air I release of COj may 
be r^id clustering of the CO 2 molecules [Eccles and Hingst, 1993], This binds the CO 2 and 
delays the 6300A aiiglow emissions that were the main observational diagnostic of the event. 
Similariy, obsNvarions of the SF6 release eiqjeriment of CRRES-at-Kwajalein suggest that the SFe 
fieezes or clusters almost totally for the initial half minute of the release, then eventually 
evi 4 )orates to effectivdy release 24 )proximatBly 10% of the intended total with a delay of almost a 
minute. Bdow we will present the evidence supporting this conclusion. 

The aiiglow measurements made by Midiael Mendillo's group at Boston University provide 
an early-time benchmark to assess the chemistry in the model. The all sky camera integrates over a 
20 second period. The measurement of the 7774A airglow cloud show no aiiglow production until 
the profile int^rated over 32 to 32 seconds Oabel 42 seconds). This suggests either (1) initial 
complete freezing (clustering) of die SF6 release or (2) complete snowplow of local plasma density 
within the release cloud for 30 seconds. The first proposed explanation is not necessarily 
uneiqiected. A theromdynamic model of the SF6 releases [Bernhardt, 1988] predicts 40% fi-eezes 
during die expansion. However, the calculation did not include partition of energy into 
translational acceleration of the gas. Much of the thermal eneigy of the SF6 gas goes into directed 
expansion out of the nozzle, thus, more than 40% may be fiozen or clustered. Soon after 
acceleration of the release gas, self-collisions will allow clustering or freezing the SFe molecules. 
Subsequent collisions with die hot ambient background will release SFe molecules, i^ch then can 
participate in attachment chemistry. Hie second proposed explanation - snowplow effects - can 
also reduce the initial 60 seconds of aiiglow production. The SFe cloud is dense in the first half 
minute (Figure 6). The density must be above 10'° cm’^ to have mean-free^ath less than 1 km. 
This situation is maintained for up to 20 s causing significant reduction of O'*' within the release. 

Even if snowplow effects can eiqilain the eariy time reduction of air^L.., Jie airglow 
measurements after 200 s indicate only 10 percoit of the SFe is available for chemistry (Figure 
11). The all sky camera unfortunately saturated at 125 Rayleighs intensity but the early and later 
time peak 7774A emissions can be cranpared to model runs with release of (a) 100%, (b) 10%, and 
(c) 1% of yapor SF«. The 10% vsqror (90% freezing) run follows die measured peak intensity after 
die camera is no longer saturated. 

The late-time chemical system was observed by a second rocket fly-through at ^proximately 
900 s after the release. The rocket trajectory and modeled ion cloud trajectory (Figure 8) were 
used to provide time history profiles of modeled IMS observatirais. The resulting picture suggests 
similar conclusions to the airglow observatirxis at earlier times. 

The resulting measured and modeled IMS profiles are compared in Figure (12). Examination 
of the three model profiles (Figure 12) suggest that 10% vapor (90% freezing) provides the best 
result fiir F', SFs’, SFj'*' doisities. This is in agreement wifo the airglow data conclusions. 
Additionally, F* as the dominant ion requires the chemical mechanism proposed by Hunton et al. 
[1987]. 

In all simulations, there were no significant amounts of SF/ produced. This differs from the 
actual experiment where SF/ was observed to have similar densities as SFe* and SF 5 *. Electron 
attachment of SF 4 is the only avenue of production included in the simulations. Very high 
temperatures are required to produce SFe* frran dissociative attachment of SFe and SFs. The 
q)paient discrq>ancy between simulation and experiment suggests that (1) there is another avenue 
of production for SF 4 *, (2) the electron attachmait rate has been underestimated, or (3) the IMS is 
measuring something other than SF 4 *. The first and second explanations rqipear less likely than 
tiien tiiird. The Mass Spectrometer used in tiie CRR£S-at-Kwajalein experiment has a resolution 
of about 3 amu. The mass of SOF 3 * (103 amu) is only 3 amu lower than SF 4 * (108 amu). SF4* 
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Peek Emission (Royteigh) 


and SOF 3 ' would be combined in the m«tsurement (D. E. Hunton, personal 

conununication). The measurement might instead predict the presence of SOF 3 *. The simulation 
results agree qualitativdy with this experimoital interpretation (Figure 13). 
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Figure 11. Peak 7774A aitglow obtained from model and All sl^ camera. Model release of a) 1.4 
xlo26 SFe molecules (no freezing), (b) 1.4xlo25 SFs molecules (90%freezing), and (c) 1.4xlo24 
SF6 molecules (99% freezing). 
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3.4 Electrodynamic Perturbation from an SF« Release 

The motion of plaana perturbations in the equatorial ionosphere requires knowledge of the 
background plasma motitm. The neutral wind blows plasma freely along the horizontal magnetic 
fields at the equator. The perpendicular to B is defined by the electric field. The bulk 
plasma moves as ExB. The global electric field system at the equator is genetaly defined by 
the interaction of the magnetic field, neutral atmosphere, and ionosphere. The laige^cale electric 
field provides the zeroth-K>rder moticm of the ionosphere. The primary ionospheric current system 
is the Solar Quiet (Sq) system. The Sq system is modeled with reascmable success using solar- 
driven, tidal-iiKxle wind models and an E rpgimi conductivity model {Richmondet al, 1976]. 

General trends of the equatorial electric fields are reviewed in Fejer [1981]. During file day 
the ionospheric plasma near the equator rises at 10 ms*^ due to the eastward electric field cieat^ 
by the Sq curroit divergences. The F region profile appears to remain stationary due to die 
ccHitinual production of plasma in the Fj region. However, after sunset the production ceases and 
the F region profile rises. The evening F] r^on chemically decays causing the bottomside 
F region ledge to steepen. The evening decay of E and F | region plasma reduces conductivity and 
causes an increase in polarization electric fields after sunset [Rishbeth, 1971]. At sunset the 
enhanconent of the downward-directed electric field drives the plasma eastward at about the same 
velocity as the thermospheric winds (100 to 200 ms*^). The F r^on horizontal drift is fiurly 
ctmstant from night to night and easily predicted fiom a knowledge of conductivities and the 
neutral wind magnitude. Unfortunately, a self-consistent electric field model is conqiuter intensive 
and does not give the desired detail needed for the evening r^on of the ionosphere. A simplified 
predictive model of the equatorial electric fields can be obtained from input of average ambient 
conditions (neutral densities, neutral winds, plasma densities, F and E region conductivities) 
[Haerendel et ai, 1991]. The calculation of large-scale electric fields (> S km) can be treated as a 
two-dimensional problem where the important ionospheric quantities are integrated along flux- 
tubes. The vertical electric field is approximately given by 





(3.8) 


where Ei is the vertical electric field, E^ is the west-to-east electric field, op is local Pedersen 
conducitivity, u„ is the zonal neutral wi^, Zff is fieldline integrated Hall conductivity, and is 
fieldline int^rat^ Pedersen conductivity [Haerendel et ai, 1991; Anderson ondMendillo, 1983]. 
Integrations are performed alotig the magnetic field line (dl) with y as a geometry factor of the flux 
tube. Figure 14 shows the integrated conductivities of the evening ionosphere. Figure IS shows 
the resulting Ei given a neutral wind of 150 ms’^, £^=1 mVm*^ and B=0.25x 10*^ T. The 
resulting horizontal motion has a shear at the F region ledge. Above the shear the velocity is 
eastward and below it is westward. 

Throughout the day the F region plasma moves upward at about 15 ms*^ Near sunset the 
plasma vertical velocity (horizontal electric field) grows in magnitude to 30-70 ms*^ for about 1 to 
2 hours. After the enhancement, the ionosphere drops rapidly (20 ms*^) and continues to have a 
downward velocity through the remainder of the ni^t. The evening enhancement of the rise 
velocity has been the focus of scientific study due to its correlation with the occurrence of ESF. 
Unfortunately, the vertical rise velocity is extr^nely variable from night to night and difficult to 
predict. 
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Figure 14. Integrated conductivities of the evening ionoq|)here. 
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Figure IS. The shear of the horizontal plasma mmion derived from the ionospheric 

conductivities shown in Figure 14. The corresponding vertical electric>field strength is shown on 
the top axis. 
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The tune histoiy of the horizontal electric field can be modeled given integrated ionospheric 
quantities and E^t) (rise velocity) [Haerendel et al, 1991]. The horizontal electric field moves the 
F r^ion iqjward. The rise alters the conductivity-wei^ted neutral wind in Equation (3.8). The 
shear feature in die hcuizootal velocity rises (and fells) with the F region. The above electric fidd 
modd conqiletes the necessary picture of the ambient conditions of neutral atmosphere, ionosphere, 
and associated electric fields. One can now introduce the chemical perturbations to the backpound 
electrodynamics and study dectrodynamics in the rc^on near the cloud. 

A two-dimensional, electrodynamic model was develc^ied to solve for the motion of die 
perturbed plasma region. The electrodynamics model is reduced to a potoitial equation with 
source terms. The source terms are determined by gravitationally-generated aiul neutral-wind- 
generated currents perpendicular to the magnetic field. The metlKxl used is similar to that 
described in Haerendel et al. [1991] and Caldr et al. [1991]. The background ionosphere is 
modded with field-line integrated quantities excluding the 30 km slice containing the main 
perturbation of the SF^ release. 

The chenustry and motion of the cloud are followed using local quantities. The reaction set 
used is sinqilified in order to reduce the number of species in the two-dimensional calculation. The 
reduced reaction set only nnodels the deep dectnm hole well. The extended chemical perturbation 
will not be modeled well with the limited chemistry. This chonical approximation is appropriate 
for determining the electrodynamic effect of tiie SF6 release. The background ionosphere is 
connected to the slab oxitaining the cloud via equipotential magnetic field-lines. 

The resulting perturbatitm of fee int^rated Pederscm conductivity shown in Figure Ifei, is 
very srruUl due to the low collisicm fiequencies at the high altitude of fee release. The largest 
electrodynamic perturbation is due to the mass loading of fee field-lines. This generates greater 
gravitational currmt in fee region of dense SFe* (Figure 166). The mass of fee SF«* plasma is still 
only a fiaction of the integrated mass of fee flux tute. However, it is sufficient to cause an 
electrodynamic perturbation. The horizontal plasma velocity perturbation is largest in the vertical 
direction (horizontal electric field perturbatimi) (Figures 16c), because the gravitationally- 
generated current is horizontal. The vertical velocity perturbation is significant up to 60 s after the 
release and should be greater in the early seconds of the release (Figure \6d). However, due to the 
rapid decay of the molecular n^ative-positive ion plasma within the cloud region, the velocity 
perturbation also will decay rapidly. This is not true of an atomic ion perturbation. 

The electrodynamics within tte SFe* cloud causes the plasma cloud to drop in altitude. A net 
drop in the plasma flux tube omtaining the ion cloud causes adjacent regions to rise. Any plume 
devek^ment comes fixnn these adjacent regitxis of the SF^' cloud. As the negative-positive plasma 
recombines, tiie net plasma loss causes fee perturbed region to slow. It could eventually rise when 
the mass in the perturbed region is less than fee surrounding r^on. 

The quick decay of the velocity field perturbation suggests that the SFe release will not 
stimulate large-scale plume structures associated wife ESF. However, the steep sides of the 
plasnut perturbation in the local release volume provide a free-energy source from which small- 
scale irregubnities can grow. The short wavelength perturbations can develop without coupling to 
the extended ionosphere. 

The horizontal velocity perturbation always remains small and the perturbation region moves 
at nearly 99% of fee ambioit horizontal velocity. The expected position of fee plasma perturbation 
can be determined by simply following the ambient motions of the ionosphere from the initial 
rdease position. 
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Figure 16. (a) The field-line integrated Pedeison conductivity, (b) The field-line integrated ion 
mass, (c) The horiaontal plasma velocity, (d) The vertical plasma velocity. The time is 60 s after 
the SF< release. 

4. DISCUSSION OF SF« RELEASES 

A thorough examination of the possible diemistry associated with an SF« release in the 
ionosphme has been performed. In making ctxnparisons of simulations and field experiments it 
was necessary to simulate (fynamics, electrodynamics and experimoital measurement procedures in 
die numerical calculatitms. Favorable comparisons between the experiment and simulation of the 
SFs release give greater confidence in interpreting the simulation and experimental data. Any 
difference s between the simulaticm results and experimental observation of the SFe release could 
reside in the oversight of significant reactions, ignorance of temperature dependencies of known 
reactions, and/or difficulties in combining complex chemistry with geometrically realistic transport 
models. The findings obtained pertaining to the release of SFe in the ionosphere as reported above 
are summarized below. 

i. The plasma perturbation moticm suggests that the lar^e-scale electric fields 
gen^ated by the perturbed field-line electrodynamic quantities are small in 
magnitude and short in duration. The perturbed plasma r^ion moves 
qrproximately with the background plasma. This suggests that the artificial 
seeding of equatorial spread F is unlikely to cause large-scale plasma plumes. 

Small-scale irregularities (< 1 km) can still form on the sharp gradients in the local 
quantities of plasma density, conductivity, etc. 

ii. The extent of the electrocfynamically perturbed region is approximately 25 km in 
radius. This is the rqgion of the deq) electron hole. The extent of the chemically 
perturbed region was determined to be very large (r > 200 km). The broad extent 
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is due to the high altitude of the release point, which allows the broad expansion of 
the SF 6 cloud. 

iii. Due to the large heat capacitance of SFt, the SFe release gas heated slowly. This 
aUows for significant levels of 7774 A airglow production, though the model 
predicted more airglow than was observed. The all sky camera does not measure 
any 7774A emission fiom the release positi(» until 42 seconds after the release 
time. This delayed emission may indicate that the SF 6 fieezes or clusters 
conq>Ietely duriug the initial e}q)ansi(Mi. The slow rise of the 7774 airglow 
sugguests an evaporation process taking place over a 1 or 2 minute period. The 
later time (~4-6 minutes) airglow measurements and the ion species profiles 
measured by the IMS (~15 minutes) best agree with about 10% - 20% of the total 
SF 6 evaporating over die initial minutes after the release. 

iv. Previously suggested SFg chemistry [Mendillo ar^ Forbes, 1982; Bernhardt, 

1984; Hunton et ai, 1987] appears to be sufficient for the initial chemistry of an 
SF 6 release (near the release point). SFe* production is large in this central regicm. 

The saturation of the all>sky camera measurement of 7774A airglow is evidence 
for the production of SFe' near the release point. 

V. P is the dominant negative ion when the cfoservaticm rocket flies through the 
perturbed r^on. The unexpected draninance of P is explained well by the charge 
interchange mechanism given in Hunton et at. [1987]. 

vi. Chemistry simulations indicated that in die late time (>5 min), oxidation of SFq is 
potentially inqiortant. Several neutral species produced later in the SFe- 
umospheric chemistry are electron attachers. SOF 3 ' is probably an important 
late-time negative ion. 

vii. SF4' densities in experimental observations are much higher than simulation 
predictions. However, the mass spectrometer resolution does not have the 
necessary resolution to separate the masses of SOF 3 ’ and SF 4 '. Therefore, the 
density of SOF 3 * predicted in the simulation may account for the SF 4 - charuiel 
IMS measuranent. 

The findiiigs of the simulation and experimoit ctmipariscais can be used to suggest directicHis 
for the continued study. The recommendations below focus on the attempt to resolve the apparent 
differences between simulation and experimental CRRES SFe release. 

While die large-scale electric fields are not perturbed gready by an equatorial release of SFe, 
small scale irr^larities can be formed and observed. Simulations of the development of the small 
scale irr^Iarities would provide insight into the timing of probable observations of the 
irr^larity onset in a positive-negative ion plasma. The dynamic and electrodynamic implications 
of die n^ativofositive ion plasma in the F region are still unexplored. The late-time modeling of 
the dqileted plasma tube can be examined in simulaticm to determine the time crmstant for refilling 
of the electron hole and plasma hole, (^estirms about the extent of the depletion along a magnetic 
field line could have possible ramifications on how big the long term electrodynamic perturbation 
is. ffithe perturbatirm has a long life, then a release of SFe cm the F region ledge prior to the pre¬ 
reversal enhancement may have a more destabilizing effect than a later release. The present 
models can suggest optimal altitudes and times for a release to seed ESF phenomenon. 

Several potentially important attachment cross sections have not been published (SFj, SOF4, 
and other molecules). Additionally, many basic chemical processes are only estimated. Products 
and rates of i(Mi-ion recombination are, in general, only estimated in this study. Associative 
deatadunoit reactions are in general unknown and left out of the present modeling. 
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Observation of the spiUia] distribution of 1304A and 13S6A emissions along with 8446A and 
7774 A emissions would provide strong indicators of the regions «^re SFe* or SFs' production is 
taking place. The 1304A and 1356A emissions must be observed fh»n a space-based platform 
such as die RAIDS satellite since the atmosphere is optically thick for these emissions. 

5. COj RELEASES IN THE F REGION 

Optical signatures of chemically induced ionospheric boles from rocket exhaust or release 
experiments have been measured and modeled for three decades [Booker, \96\,Mertdillo and 
Forbes, \91%-,Zinn and Sutherland, I9i0; Mendillo and Baumgardner, 1982]. Releasing CO 2 in 
the F r^on allows for observations of CO} chemistry in a time-dependent regime which is difficult 
to study in the laboratory. The Red Air I (Release Experiments to Derive Airglow Inducing 
Reactions) Program released 9 kg of CO 2 at 2S3 km and again at 379 km altitude in the mid¬ 
latitude nighttime ionosphere from a single rocket. The strategic release altitudes help to isolate 
key chemical processes to determine reaction rates of COa/ionospheric chemistry. 

The presence of CO 2 in the ionospheres of Venus and Mars as well as the exhaust effluents of 
some rocto engines make the study of CO 2 reactitms an important endeavor. CO 2 is an important 
perturbing chemical in the high ionosphere (z> 200km) of tte earth. It is not naturally present in 
any significant density but the ambient plasma is perturbed strongly by the release of CO 2 . This is 
because the reaction rate of O'** with CO 2 is 2 atul 3 orders of magnitude greater than O'*' with 0} 
and N}, respectively. A CO} release in the F region ionosphere rapidly changes the plasma friom an 
atcnnic to a molecular ion plasma, thus increasing die rate of ion-electron recombination. 

The primary reactions of CO} in the F region are. 

COj+0*->0;+C0 (5.1) 

COj+0* ->00^+0 (5.2) 

The total rate of the combined reactions (5.1-5.2) is a{^roximately 9.4x 10*^*^ cm^ s'^ and does not 
dqiend greatly on temperature. However, the branchir^ of the C0}-0'*' reaction does vary strongly 
with leactimi energy [Fehsenfeld et al, 1966; Paulson et al, 1966; Dunkin et al, 1968; 
Mosesman and Huntress, \910-,Johns€netal., \910-, Lindinger et al, \91A-,Jqffe and Klein, 
1974; Viggiano et al, 1990a and 19906]. Implications of the branching energy dependence for 
ionospheric release of CO 2 have been discuss^ in Caledonia et al [1987] and Hunton et al 
[1990]. At low energies, <1 eV, the charge interchange branch is dominant. In rocket release 
experiments, this C0}-0'*' interchange produces O}'*' and, eventually, 6300A airglow. At higher 
energies the charge exchange branch dominates. The same branching trends are observed with 
both vibrational energy and translational energy (A.A. Viggiano, private communicatimi). 
Additionally, the pr^ulation distribution of the vibrational states of the produced O}'*' varies with 
tonperature [Viggiano et al, 1990a]. This results in a change in the vibraticxial population of O}'*' 
and altered branching ratios of dissociative recombination [Guberman, 1988]. 

Due to the rarefied conditions of an upper atmospheric release of CO}, NLTE (non-local 
thermal equilibrium) conditions hold. Under such conditions the vibrational state populations of 
molecular neutrals and ions can depart fi’om Boitzmarui distributions. Chemical rate coefficients 
obtained under LTE conditions may be a misrepresentation within moctels of chemical release. 
Explicit vibrational dependencies in rate coefficients and products are potentially important in 
modeling such NLTE environments. 
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Tlie molecular ion plasma produced in a CO 2 release has observable effects. A large release 
might lead to a long-lived plasma depletion through increased recombination rates and to a 6300A 
aiiglow enhancement through the radiative decay of 0 ( 1 D) produced in the rearnibination of O^'^: 

Oj+e - >0CD) + 0 (5.3) 

and possibly in the recombination of CO 3 '*' 

CO^ +e- >0(^0} *-CO (5.4) 

The 6300A emission perturbation and the electron density depletion provide signatures of the CO 7 - 
ionosj^ric diemistry and tirtte-dependent airglow chemistry. 

COj-ionospheric chemistry has been used to model 6300A airglow production in C 02 releases 
firom exhaust plumes [Mendillo and Baumgardner^ 1982]. There were indications that the 
observed 6300A intensity was less than the value predicted by the chemistry model. 

The Red Air I Campaign attempted to place CC )2 releases above and below the F peak. For 
these two release points, the electron densities wmild be similar but the neutral background would 
give different characteristics in diffusion and chemical constituent concentratiens. Measurements 
of the time-dependent airglow perturbations at two separate altitudes in the F region have the 
potential to provide definitive information on unresolved questions of nighttime airglow chemistry 
and CO 2 diemistry. If published C 02 - 0 '^ reaction rates are accurate, then CO 2 releases at several 
altitudes can be used to find quenching rates of 0(^D) by ground state O. At altitudes below 250 
km N 2 quenching quenching of 0(^D) dominates. Above 290 km radiative quenching dmninates. 
Near 270 km altitude O quenching of 0(^Z)) may be important. This quenching reaction, first 
suggested by Yee andDalgamo [198S], could alter pre^cted airglow profiles for daytime, 
ni^ttime and auroral ionospheres. Some investigators suggest that quenching by 0 is not needed 
to e?q)lain experimental observations [Link and Cogger, 1988] (also see. Bates [1989]). 

This investigation examines the chemistry and thermodynamics of a cold CO 2 release in the F 
region of the ionosphere. The direct comparison of model results and Bostmi University all sky 
camera observations made during Red Air I erqrerinKnt. The discrepancy between model results 
and observations of the Red Air I releases indicates that either (1) CO 2 -O'*' reaction rates in the 
literature are much too first, or ( 2 ) complex diermodynamic processes are occurring that degrade 
the intended sirrqrlicity of the experi ^lent. Several questions are examined in numerical simulations 
of the Red Air I release to explain the observed airglow signatures. First, what processes can 
produce the less-than-expected airglow emission fi-om the €(>2 release region? Second, what 
processes can lead to a binuxlai decay curve observed in the high altitude release? 

Because the primary observable in the Red Air I experiment is 6300A emission fi-om radiative 
decay of 0(^D), it was necessary to examine nigh^low chemistry models. The reaction set 
necessary fin- ambient 0(^D) productimi and loss is reviewed elsewhere [Torr and Torr, 1982; 

Link and Cogger, l9iS; Solomon et al., l9iS; and Solomon and Abreu; 1989]. The backgroimd 
6300A airglow during Red Air I was measured and is modeled as an additional baseline input into 
the study of CO 2 ion-chemistry. Simulations use two separate airglow chemistry sets — airglow 
diemistry suggested by Link and Cogger [1988] and the chemistry suggested by Solomon and 
Abreu [1989]. In the near steady-state conditions of the ambient nighttime airglow, ground 
measurements of the integrated airglow intensities caruiot distinguish between the two chonistry 
sets [Link and Cogger, 1988]. However, the time-dependent, altitude-specific effects of the COi 
releases forces one to consider the differences between the assumed rates. 

SecticMi 2 is a description of the Red Air I experiment. Section 3 of the paper discusses 
ambient airglow chemistry and implications of the Red Air I optical observations r^rted by 
Michael Mendillo [perstmal communication]. Section 4 examines the results from 
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muldduiiensional release modeling. Section 5 presents the simulations of the Red Air ICO 2 
releases using several more complex chemical models. Conclusions of the p^^>er are in Section 6. 

6 . RED AIR I EXPERIMENT DESCRIPTION 

The Red Air I Can^gn attenqited to place CO 2 release above and below the F peak. For 
these two release points, the electron densities would be similar but the neutral background would 
give dififeieot diaracteristics in diffusion and chemical constituent axicentrations. A stuile rocket 
carried both release canisters. One was released at 2S3 km altitude (2:29 ITT) and the other at 379 
km altitude (2:32 OT). Table 6 provides the geosolar conditions for the night of the Red Air I 
release ejqieriment. 


Table 6. Release Conditions. 


Launch 
10.7cm Flux 
< 10.7cm Flux> 
Sunspot Number (Ra) 
Geo-Longitude 
Geo-Latitude 
^Index 
Kp Index 
Release Altitudes 
Local Solar Tiffle 


3 April 1989 
184 
196 
218 

284.6® E 
37.8° N 
21 
27 

379 km and 253 km 
21:30 


The Millstone Hill radar measured the electron density profile, electron temperature, and ion 
ten^ierature during the time of the airglow measurmnents. The MSIS90 neutral atmosphere model 
[Hedin, 1991] predicts an exospheric temperature of 1200K for the gec^hysical parameters of the 
Red Air I releases. The neutral exospheric temperature of the MSIS model atmosphere was set 
equal to tiie measured itm tonperature of the R^ Air I experiment (1030K) [Mendillo, personal 
communicatitxi, 1991]. 

All sky camera observatitms of the 6300A airglow enhancement showed an enhancement of 
75 Rayleighs over the background for the upper release of CO 2 . The enhancement was immediate 
(< 1 minute) and decreased rqiproxinudely as a simple decay curve, though the decay curve had a 
bimodal characteristic suggesting the presoice of other processes b^ond simple C02-ionospheric 
chemistry (Figure 17). 

The observation of the low-altitude release showed a peak enhancement of 25 Rayleighs over 
the background. The enhancement time history of the lower release had a more pronounced hump 
diaracteristic. The double rise of airglow at the low-altitude release positions also suggests that 
processes odier than simple CO 2 chemistry may be involved in the perturbation evolution. Note 
that the sec(»d pronounced rise in the low-altitude release measurement corresponds to the release 
time of die high-altitude release. 

It is important to examine release characteristics, species transport, and chemistry in order to 
fully interpret the Red Air I experiment. Rather than being a definitive experiment to quantify 
questions cm ambient airglow and CO 2 ion chemistry, the experiment is instead an investigation 
into unexpected thermodynamic and chemical processes. 
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Figure 17. Aiiglow measurements from a 41^ 4 degree region centered on the upper release 
(triangles), centered on the lower release (squares), and on the undisturbed background (circles). 
The badcground measurements are smoothed. 


7. NIGHTTIME BACKGROUND AIRGLOW 

While metastable species are only trace crmstituents in the atmosphere, they have significant 
roles in the diemistiy and diermal budget of the atmosphere. Lrnig radiative lifetimes permit 
conversion of electronic energy to kinetic or vibraticmal energy through quenching processes. In 
addition, the Icmg-lived metas^le species allow for energy transport away fhxn the r^on of the 
metastable species firrmation. The metastable species 0(^D) provides one of the strongest features 
in aitglow and auroral spectra. Due to the relative simplicity of the chemistry of the nighttime mid¬ 
latitude icmosphere, nightglow measurements from this region can help resolve rate coefficient 
discrqrancies of the 0(^D) sources and sinks. The ambient reactions used to model airglow 
production are given in Tables 7 and 8. 

While Hays et al. [1978] used Atmosphere Explorer (AE) satellite nightglow measuronents to 
determine quenching rates and yields of 0(^D), Cogger et al. [1980] used the chemistry of Hays et 
al. [1979] to predict nighttime 6300A airglow and compared the prediction to ground-based 
measurements. They concluded that the chemistry predicted airglow intensity at twice the observed 
intensities fixxn 0(^Z)) decay. Through the inclusion of several new rate determinations. Link and 
Cogger [1988] reexamined the ground-based woric of Cogger et al. [1980] and resolved the 
discrepancies. However, some issues of 0(^D) nightglow chemistry remain (for a fell discussion 
see Link and Cogger [1988]). 

The primary issue highlighted in Link and Cogger [1988] (henceforth LC) is the quenching of 
0(^D) by ground state atomic oxygen. The nightglow chemistry put forward by Solomon and 
Abreu [1989] (henceforth SA) includes a non-zero rate for this quenching reaction. The 
differences between the IC and SA chemistry sets center on the use of 6.81 x 10'^ s*^ and 9.41 x 10 * 










3 s'^ reflectively, for the Einstein coefiBcient of the 0(^D)-0{^F) transition. To agree with 
observed aiiglow values SA included the quenching of 0 (^D) by atomic oxygen 

0('X>) + 0(*i»)- >OCP) + Oi^P) (7.9) 

with a suggested rate of 8x10*^^ cm'^s*^ LC stated that the quenching of 0(^Z)) by radiative 
decay (3.9) is not needed if a smaller Einstein coefficient for 0(^£>) radiative decay is used. 

A steady-state model of the ionosphere was used to investigate the background int^rated 
aiiglow associated with the Red Air I release experiments. The model assumes the ionospheric 
constituents are in chemical equilibrium. The ion constituents can be determined from a measured 
electron density profile, given neutral density profiles, and a set of rate coefficients for the 
reactions listed in Tables 7 and 8 . 


Table 7. Nighttime airglow chemistry. 


DISSOCIATTYE RECOMBINATION 


Nj +e—!^N + N 

At, =1.9x10-’( 7 ;/300)““‘' 

NO* +e—i^N + 0 

=4.0xl0-’(r,/300)““’ 

Oj +e ‘‘*-40+0 

itj=2.0xl0-’(r,/300)^’ 

CHARGE EXCHANGE 


N* +0—+N 

= 2.0 X 10-*’ 

N* +NO-^NO* +N 

*5 = 6.0 X io-'®(r„/3oo)“” 

N*+Oj—J^O; +N 

= 3.0 X 10-'° 

0 * +n(’Z))—!^N* +0 

= 13 X 10-'° 

0 * +NO-^NO* +0 

it, =8.0x10-'’(r„/300)“®^ 

o*+Oj —+0 

it, =2.0xlO-"(7;/300)'®’ 

Nj +0—^0* +N 2 

it,o = 1.0x10-" ( 7 ;/300)'®" 

N; +N0—!^N0* +Nj 

it,, =3.3x10-'° 

N*+Oj^^O*+Nj 

it,j =5.0xl0-"(r„/300)'^* 

0 ^+n(*o)—^N *+02 

it,, = 23 X 10-'° 

CHARGE INTERCHANGE 


N* +0j—^0* +N0 

it,j =4.0x10-" 

N* +0a ->N0* +0 

Jt,g = 4.4 X 10-'’ exp[-3220 / T,] 

0* +N; ->N0* +N 

it,7 =12xl0-'’(7;/300)''“ 

n; +0—^NO* +N 

= 1.4 X 10-'°(r„/300)'”^ 

0 ; +N—^NO* +0 

it,, = 1.8 X 10-'° 
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Table 7. (continued) Nighttime airglow chemistiy. 


N+NO—+0 

kjo= 3.6x10’" 

N + Oj—!^N0 + 0 

kj, = 4.4 X10’’^ exp[-3220/ T,] 

N(’D) + NO—!^Nj +0 

kjj =6.0x10’" 

N(*£>) + 02—!^N0 + 0 

kj,=8.8xl0’‘*(7;/300)®^ 

N(*D) + 0—!^N + 0 

k„ = 1.9 X 10-'^ 

N(*£>) + Nj—!^N + Nj 

kjj = 9.4 X lO"'* exp[S 10 / T„ ] 

N(*D) + Oj—!^N + Oj 

kjs =53x10’’^ 

N(*Z)) + e—!^N + e 

kj7 =3.9xl0’'°(r./300)°*' 

N(*P) + 0— J^»->N + 0 

kjg =13x10’" 

N(*P)+e—!^»->N+e 

itj9=2.8xl0’®(r./300)°* 

0(‘D) + 0—!^0 + 0 

it3o= 0.0-8.0x10’" 

CK'/5) + Nj—^0 + Nj 

k„ = 2.0 X 10’" exp[l 11 . 8 / 7 ;] 

0('D) + 0j— + 

k,j =2.9xl0-"exp[673/7;] 

0('£>) + e— ^0 + e 

kj, =2.9x10’'® ( 7 ;/300)°®' 


Table 8. Radiative quenching. 


N(’D)- 

-2»-^N + hv(5200A) 

a, = 1.07xl0*^/5 

N('P)- 

-2i->N + hv(3455A) 

02 =5.4x10~^/5 

N(*P)- 

-^N(^jD) + hv(10400A) 

a3 = 7.9xl0“^/j 

0(D)- 

0 + hv(6300A, 6364A) 

04 =6.8-9.4x10 "^/5 

OCS)- 

-2i-^0 + hv(2972A) 

a5 = 4.5xl0"^/s 

0('i)- 

-^0(!D) + hv(5577A) 

Ofe = 1.06/j 


The background neutral atmosphere is detennined by MSIS90 [Hedin, 1991], which provides 
die altitude profiles of N, 0, N 2 , NO, O 2 . Table 6 lists the geo-solar conditions of the Red Air I 
experimoit. The electron density was obtained from Millstone Hill radar measurements (M. 
Mendillo, private communication). The main ion species NO'*', O’**, O 2 ''' are determined through 
steacfy state assumptitxis that lead to relationships similar to that in Link and Cogger [1988]. 
Finally, minor species and metastable species N(^D), 0('Z7) and 0('S), and N'*' are determined 
firom the neutral and ion species mentioned above, again using steady state assumptions. 

The observed background 6300A airglow for the Red Air I experiment was detennined from 
uiqiertutbed portions of the all sky camera images. The measured background airglow was near 
90 Raylei^ fiir a zenith look. This is equivalent to 108 Rayleighs with an atmospheric extinction 
correction of 1.2 [Baumgardner, personal communication, 1991]. 
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The model deteiminatiiMi of integrated 6300A aiiglow intensities were goierally much higher 
than the measured intensities. The aiiglow chemistry of SA and LC predicted aiiglow intensities 
that were approximately two times the observed intensities. There are several parameters of the 
mod^' can, if altered, reduce the modeled aiiglow, e.g., error in diemical rates, exosphoic 
tenij are input into MSIS90, error in the measured electron density profile, etc. The lower 
atmospheric tenqrerature of 1030K was used for determiniitg the atmosphere density profiles and 
reaction rates. This alone reduced the predicted aiiglow by 50%. To match the measured aiiglow 
value, the rate of the reaction 

0*+0, - >01+0 (7.10) 

was assumed to be lower than the rate given by Albritton et al. [1977]. The experimental 
conditions of Albritton et al. [ 1977] de-coupled the kinetic motion energy of the reactants fiom the 
O 2 ''' vibrational energy. Chen et al. [1978] measured the rate of reaction (7. 10 ) wifii true thermal 
distributions for the reactants and O 2 ''' internal energy. This resulted in lower reaction rates than 
Albritton et al. [1977]. Unfortunately, the Chen et al. [1978] measurements were fixxn 300K to 
900K. The observed aiiglow can be reproduced in the model with no alterations of other rate 
coefficients given in the SA chemistry set when the formula suggested by Chen et al. [1978] for 
reaction (7.10) is used with a temperature of 1030K. 

For the electron density profile and plasma temperature measured for the Red Air I 
experiment, the steady*state chemistry model predicts 6300A aiiglow values of 110 Rayleighs and 
140 Rayleighs for SA and LC, respectively, oxnpared to the measured 6300A airglow value of 
108 Rayleighs. SA provided excellent agreon^ with the experimental measurement. The 
production rates of 0(^£>) are the same for both chemistry sets (Figure 18). LC assumed no 
quenching of 0 (^f>) by 0 and the slower radiative decay for 0 (^Z)). 


0(1 D) Production Rotes 



Production Rote (cm“^ s~') 


Figure 18. The production profiles of 0(^D) for both LC chemistry set and SA chemistry set. 









Figure 19a shows the loss rates of 0(^f>) for LC chemistry. A OOD) quenching rate of 4x 
10'^2 cm*^s*^ was needed to bring LC chemistry into agreement with the Red Air I badcground 
aiigiow measurement. The loss rates of 0(^f>) by various mechanisms are plotted in Figure 196. 
for SA chemistry. It is seen that 0(^£>) quenching by O has only a narrow altitude range of 
domination as the main loss mechanism, but it has a broad altitude range of relative inqwrtance. 
The first low altitude release occurs .fortuitously, precisely at the altitude of the peak O quenching 
of However, d^rading effects reduce foe reliability of the lower release of Red Air I for 

quoiching rate determinations. These degrading effects are examined below. 

Figure 20 compares the normalized int^rated airglow as measured by a zenith loddng rock^ 
instrument for SA and LC. There is a measurable difference if foe observing rocket reaches an 
altitude of 300 km. The present ambiguity of tte ground based measuronents of ambioit 6300A 
airglow leaves no indication as to the preferred Einstein coefficient of O^D) decay and the 0(^D) 
quenching rate by O. However, foe time-dependent emission perturbation induced by the Red Air I 
CO 2 releases help resolve the quenching and radiative issues of 0(^£>) chemistry. 



Figure 19. The loss rates of 0(^D) by various mechanisms using (a) LC chemistry set and (b) 
SA chemistry s^. 
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Normalized Airglow 

Figure 20. The normalized integrated 0(^D) airglow as measured by a zenith looking rodcet 
instrument for the LC chemistry set and the SA chemistry set. 


8 . TWO-DIMENSIONAL SIMULATIONS OF RED AIR I 

The CO 2 releases at 2S3km and 3791an gave 6300A airglow enhancements that differ with 
each other in magnitude, rise time, and decay shape. The differences are not simply explained by 
the chennistry results in the above section. Modeling the Red Air I experiment requires simulation 
of die dimnical and dynamic effects of the release and the observation methods of the airglow 
measurement. As will be seen, the getmietry of the line-of-sight integratitm of the airglow played a 
role in the rruignitude and time history of the 6300A measurmnent. 

The Red Air I releases of CO 2 into the mid-latitude nighttime itmosphere are modeled with a 
two-dimensicmal finite-difference model. Cylindrical geometry, with the axis of symmetry being 
altitude, was used to reduce the three-dimensicxial transport problon to a two-dimensional 
numerical computation. It is assumed that the magnetic field is vertical and the horizontal rocket 
motion is negligible. The diffesion limit is used for species transport. The ambient species for the 
initialization of the background within the two-dimensional model are provided by the steady-state 
model discussed in Section 7. C02 is introduced at the release altitude within the mocfel and the 
diemistry and diffesitm of the ambient and caib<m oxide species (C02, CO 2 '’', CO) are followed in 
time. 

An average Rayleigh value is obtained fi'om simulation data to compare with Red Air I data. 
The simulated airglow peak is rqrproximately four times too large for both the SA and LC reaction 
set (Figure 21). This suggests that a chemical or thermodynamic process is missing or incorrectly 
represented in the model. In addition, the shape of the enhanced airglow curve of the simulation is 
a simple decay structure and not the double hump of the Red Air I data. 

Figure 22 displays the simulated lower release airglow. The magnitude and the sluqre are not 
modeled well in the lower release. The two releases of the Red Air I experiment were performed 
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from a smgle rocket. This allows for the possibility of mixing the aiiglow signatures of the two 
releases spatially. This mixiog is shown below to be important only for the lower release. It is 
interesting to note that frie SA chemistry produced the brightest aiiglow in die upper release and the 
weakest response in the lower release. Tbis is due to the peak quenching of 0(^^) by O at the 233 
km altitude release. 



Figure 21. The modeled and measured 6300Aairglow from the center of the ujqier release. Tlie 
modeled airglow is obtained from the SA chemistry set (circles) and the LC ch^stry set 
(squares). The LC set airglow curve is divided by four to show the approximate match with the 
data (diamonds). The Red Air I data frpom the iqiper release is plotted with triangles. 
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Figure 22. The modeled and measured 6300A airglow from the center ofthe lower release. The 
modeled airglow is obtained from the SA chemistry set (circles) and the LC chemistry set 
(squares). The LC set airglow curve is divided by four to show the approximate match with the 
data (diamonds). The Red Air I data from the lower release is plotted with triangles. 
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To detefmine if the aiiglow perturbations of the separate releases were "mixed", the airglow 
histories of tibe two simulations are linearly a dded using the observational geometry of the all sky 
camera. The rdeases were sqparated enough to make the linear adding of the airglow effects 
reasonable. The CO 2 of the lower release may have expanded up to the higher altitude of the 
second release and vice versa. Figure 23 shows die ad^ effects of the two releases at die airglow 
measurement point of the high and low altitude releases. The eiqiaiisioo of the lower release cloud 
is st(^q;ied abrupdy at the lower altitude. The lower release does not affect die upp^ release 
r^ion. The shape of the tune decay of die airglow enhancement measured at die upper release 
point cannot be explained by the mixing ofdie simulations ofthe two releases. Aseparate 
chemical process must be at work in die upper rriease. This is examined in the next sectimi. 

The hitler altitude release has a definite influence on the measurement of the lower release 
eiqieriment. The upper release expands and fidls into the r^on of the lower release point. The 
IM Air I airglow measurement has a plateau near the fiont of the time-airglow curve. This could 
be the mixing r^on of die upper and lower cloud effects. This is unfortunate since this is the 
r^on vhere O quenching of 0(^D) is important. The lower altitude intensity predicticms with LC 
and SA chemistry are both too high by a fiictor that is larger than the upper release predictions. 

The ratio of the peak intensity of the upper to the lower releases is approximately 3.8. The LC set 
produces a ratio of 2.0. The SA set with O on 0(^£>) quenching pnx^ces a ratio of 2.S, which is 
closer to the measured value. The absolute intensities are affected by fireezing of CO 2 up<m 
release. Perhaps the upper and lowm releases have different release characteristics duU allow 
different percentages of CO 2 to participate in the airglow chemistry. This would invalidate any 
conclusion based on the ratio of upper to lower airglow intensities, but the ratios indicate that O 
quendiing may be present. Two remaining questions are addressed in the fiiUowmg section; (1) 
why does the model over-predict the airglow, and (2) what produces the hump in the upper release 
airglow signature? 



Figure 23. The measured airglow intensities for both upper (empty triangles) and lower releases 
(solid triangle) are compared to the predictions of the S A chemistry (dashed line) and the LC 
chemistry (solid line). The airglow predictions have been divided by a constant so that the upper 
release measured and modeled airglow intensities are the same. 
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9. COj CHEMICAL INVESTIGATIONS 

Conqnitatioiis were perfomied using a sin^-cell Lagrangian chemical-kinetics code to 
investigate die possible differences in the chemistry sets of standard CO 2 reactions, explicit 
vibrational C02-02'^ reactions, clustering COi reactions, and the freezing and subsequent 
sublimation of C02 reactions. The removal of CO 2 by diffusion is also included in tlx model. 

The focus of this investigation was to determine chemical processes that may produce the 
hump observed in the hi^ altitude release 6300A emission profile. It should be that without 
spatial and dynamic modeling, it is difficult to (fotain absolute comparisons with airglow 
measurements. 

The initialization of model parameters was based upon the actual release conditions. Each 
release consisted of 9 kg (1.24>(10^^ molecules) of CO}. Based on data from the all sky camera, 
the release at high altitude is assumed to have quickly expanded to a cloud having a rachus of 2S 
km. At this size, assuming homogeneity, the cloud CO 2 density was estimated to be 10° molecules 
per cubic cratimeter. The rate that CO 2 left the system via diffusion in the model was estinufted to 
be 10‘^ cm^s‘1. 

9.1 Standard CO 3 Chemistry 

hi tile first attempt to itKxlel the release, the standard CO 2 reaction set outlined in Section 5 
was used. Fimies 24 a and A rqiresent the time histories of CO 2 and 0(^Z>) density. From the 
density of 0(*D) a rough estimate of tiie peak oihancement over the background can be made. 
Assurning a 50 km rcgimi of relatively dense CO 2 , tiie peak enhanconent from this model was 
estimated to be 148 Rayleighs. This is twice the diserved peak. From the 0(^D) profile shape it 
is clear tiiat tiiis simple model does not produce the hump in the observed emissions. It was 
determined tiiat tiie standard chemistry dissipates CO 2 too rapidly to create the hump. This 
suggests that there are other, previously unsuspected processes that affect the CO 2 chemistry. 

9.2 Explicit Vibrational CO 2 Chemistry 

The inclusion of eiqilicit vibrational chemistry was an attempt to model the effect of heating 
COi vibrational levels on its decay. Introduction of explicit vibrational states of CO} (A. A. 
Viggiano, private communication) and O 2 '*' [Guberman, 1988; Queffelec et al., 1989] into the 
chemistry is examined to allow for non-local tiiermal equilibrium effects on the simulated 6300A 


airglow. Reactions (1.1) and (1.2) become 

CO 2 (0000) + 0'"- > Oj + CO 

(9.1a) 

CO2(0110)+O'"- 

-^Oj+CO 

(9.1b) 

C02(0000) + 0^- 

—^COj+0 

(9.2a) 

CO2(0110)+O"^- 

—>coJ+o 

(9.2b) 


The rates of reactions (5.la,6) and (5.2a,6) are given as a function of vibrational states of CO 2 
using uiqmblished data of A.A. Viggiano (private communication). There has been no attempt to 
extrapolate the data up to the high ionospheric tonperatures, thus, data for T^^q^ = SOOK is used. 
Rates for vibrational quenching and excitation of CO 2 are taken from Archer [1990]. 
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Brandung of the COrO'*' reaction path varies with vibrational level. Ground state CO 2 
producing O 2 '*' (9- has a branching ratio of 0.96, leaving reaction path {9.2a) relatively 
unimportant to the system. COjCOllO) producing O 2 '*' (9.1b), however, controls only 87% of the 
reactioiis, making the producdoo of CO]'*’ (9.2b) in^rtant thou^ still small. In ad^on to diese 
changes, all low enogy reactions (E < 0. leV) involving CO] were updated to include these explicit 
levds. The reactions wme taken from .drc/ier (7990). 

The production rates of vibrational states of O 2 ''' are taken from Viggiano et al. [1990a]. 
These rates are based on data taken in the temperature range of 88K to 300K. This is a limitidioa 
on die NLTE investigation but die results will indicate trends, (^endiing rates of 02'*'(m) are 
estimated. Further labmatoiy studies are necessary to provide proper rates for die dtemioal 
models. 

Figures 2S a and b represent the time histories of important species during the model 
simulation. It is clearly seen that 0(^Z)) production is not changed in any significant way fixxn the 
standard CO] model results. The decay slu^)e and time constant are virtually identical to that 
previous result. The plot of CO] doisity shows its decay is also unaffected. The peak 
enhancement over the background is estimated to be slighUy more than a foctor of two too great 
and the characteristic hump has not been created. This implies that the non-equilibrium effects of 
vibrational excitation of CO] are not responsible for the secondary peak of 0(^0) production diat 
produces the observed emissions hump. 

Conqmtations using the two diffbrent chemistry sets of CO] (with and without vibratiaoal 
states) give the same results. There is presendy no advantage to using oqilich vibrational levels of 
CO] and O]*^ in the diemistry. This does not say that vibrational level diemistiy is not important 
in the modeling ofthe NLTE conditions in CO] releases. A possible oqilanation is foat die rates 
for reactkms omtainii^ CO] and O]*** were obtained fixim low temperature measurements 
(T<545K and 300K, respectively). Laboratory eiqrerimaits with higher temperature environments 
could give reaction rates that would alter results of subsequent simulations. 

The results indicate duU too much 6300A aiiglow is predicted by the best simulation effort. 
Possible reasons for the low 0(^73) production in ^ Red Air I release are: 1) COj freezing or 
clustering at release, 2) error in the MSIS86 neutral atmosphere, 3) error in the measured electixm 
density, 4) slower charge interchange rate of CO] with O'**. 

9.3 CO3 Cluster Chemistry 

The diemistry kinetics code was used to investigate the impact of CO] clustering (or freezing) 
and subsequent evaporation (or sublimation) on the induced 6300A aiiglow produced by a release 
of CO] in die ionosphere. Two eiqierimental observations in the Red Air I diglow nreasurements 
suggest that clustering or freezing plays a role in the CO] release chemistry. First, the observed 
6300A emissions are several times less dian predicted emissions. The reduction of modeled 6300A 
aiiglow magnitude can be achieved if CO] fieezes or clusters. Simple ideal gas calculations 
pmlict that 7% of the CO] would freeze in an ionospheric release and be lost to the experiment 
[Berrthardt, 1987]. Because of a huge dipole moment, the CO] molecule is for from ideal. Van 
der Waals fiirces should be taken into account in the freezing calculatimi including Joule- 
Thonqison cooling and kinetic energy loss to directed translational acceleraticm. Thermodynamics 
calculations suggest dud more than 50% of the total CO] release may cluster or freeze. Tte 
proposed productirm of clusters was corroborated unanimously by several clustering experts. 
Indeed, huge clusters (n > 100) were likely (Olaf Echt, University of New Hanqishire, personal 
communication). This would reduce the amount of free COj. 
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MODEL n. VIBRATIONAL CHEMISTRY 



Figure 25. (a) Time histoiy of CO 2 vibrational level density for high altitude release, explicit CO 2 
vibrational state chemistiy set. (b) 0('D) density time histoiy for high altitude release, explicit 
CX)] vibrational state chemistiy set. 
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The second reason for considering freezing and clustering is the presence of a hump m the 
decay of the 6300A aiiglow signature. The hump might imply the slow release of free CO 2 from 
the frozen 01 clustered CO 2 . 

It is unclear how die clusters actually decay in collisicm, requiring that two cases be 
investigated. The first case assumes that after sufficioit energy has been imparted to the cluster 
via cdlisions, a single COj molecule breaks off from the larger group. The seccmd case models the 
ejecdon of a smaller cluster from a larger heated group. 

For both of the models it is assumed that cluster breakup occurs only due to collisions with 
relatively pc^ulous atomic oxygen and diatomic nitrogen. An estimate of the reaction rates for 
these processes is made by taking into account die colliding particle cross-sections and the oieigy 
of the colliding particles. The basic equation to determine die reaction rate is then 




(93) 


where oq is the collision cross section, k is Boltzmann's constant, T is temperature, and M is the 
particle mass. The subscript T refers to values for the CO 2 clusters and “2" represents values for 
either O or N 2 . For this calculation, the temperature of the cold CO 2 is assumed to be 250K and 
the hotter ambient species ate 1030K. The masses of the individual particles are 7.304x 10'^^ kg 
for CO 2 ,4.648x lO*^^ kg for N 2 , and 2.6S6x 10'” 1^ for 0. The collision cross section is assumed 
to be circular and is determined by 

(To^^(n+r2f (9.4) 

where r is the particle radius. The radius of CO 2 is 2.2SA, the radius of N 2 is 1.9A, and the radius 
of 0 is I. 3 A. The efiective radius of the CO 2 clusters, rj, varies with cluster size, n, in the 
following ituuiner 

n = «*^C02 (9 5) 


where is the effective radius of a single carbon dioxide molecule. 

It is assumed that a pardcle colliding with a CO 2 cluster transfers energy only to a single CO 2 
molecule in that cluster, and this energy is then transferred to the rest of the cluster via weak 
cluster bonds. Whoi enough energy is absorbed into the cluster then a particle will “boil” off 
This energy transfer collision frequency is 


= 


f M, 


\My+M2J 


U 


(9.6) 


^lere the mass is no Irniger the cluster mass but the mass of one C02 molecule. The mass 
weighting in equation (9.6) for energy transfer collisions is 0. IS for N 2 collisions and 0.07 for O 
collisirms. 

It is also assumed that the binding energies of the clusters are small and a single collision will 
remove molecules from the cluster. The rates for reactions developed with this model are found in 
Figure 26 below as a functitm of cluster size up to 10. 


44 







CLUSTER COLLSION REACTION RATES 



FigufC 26. Estimated reaction rates for breakup of clusters due to collisional heating. 

9.3.1. Single Particle Removal 

The diemistty for this model is identical to the standard CO 2 reaction set with the addititm of 
reactions describing the breakup of the CO 2 clusters by ejection of a single CO 2 molecule. These 
reactitns are of the form 

(CO2 )n + M->(C 02 )„_i + CO2 + M ^97 

vdiere M is either O or N 2 . The rates for these reactions are determined as described above. 

Figures 27 a and b display the time histories of 0(^D) and CO 2 density for this model. Fifty 
percent of the total CO 2 cloud density is assumed to be clustered upon release as (C02)io clusters. 
Larger cluster sizes (n=30,60,90) were also used but n=10 produced results that better reflected 
the 6300A decay curve. Figure 27a shows that the simple breakup cascade of CO 2 clusters down 
to single molecules does, in foct, delay the formation and relaxation of 0(1/)). From that resulting 
0(^D) history a peak enhancement of 92 Rayleighs is estimated with no hump predicted. The 
predicted intensity is now very close to the observed 75 Rayleighs. It is clear t^ this model does 
not delay the introducticni of CO 2 into the system nearly long enough to create a secondary 0(^D) 
density peak wdiich would lead to the obsen^ doisity curve. However, the reduced peak 
emissions foom relaxing 0(^Z7) indicates that a significant percentage of the released CO 2 may be 
clustered, verifying a critical assumption made. If the second assumption of small cluster sizes is 
valid, the siu 4 >e of the oxygen curve suggests that boiling off one CO 2 molecule at a time finom a 
heated cluster is not the correct decay mechanism for the process. The supply of a significant 
amount of CO 2 to the system must come after the initial release and it must come fiisf enough to 
make a secrmdary peak in the carbon dioxide decay profile. 
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MODEL ma. CLUSTER CHEMISTRY 
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Figure 27. (a) CO 2 cluster density time history for high altitude release, single particle removal 
model: (C02)o + M (C02)n.i + CO 2 + M. 50% of total CO 2 released as (C02)io- (b) 0(^D) 
density time history for high altitude release, single particle removal model; (€02)^ + M 
(C02)n.i ^ *°bil CO 2 released as (C02)io 
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9.3.2. Small Cluster Removal 

In an attempt to more accurately model the decay of CO 2 clusters, reactions of the fidlowing 
finm were added to die standard CO 2 reaction set 

(C02)„+M-►(C02)„_2+(C02)2+M (9g) 

where M again rqiresents O or N 2 . These reactions have the same reaction rates as their 
counterparts in the single particle removal model. CO 2 is then released via 

(€02)2+M->COj+COj+M ( 9 . 9 ) 

Again it was assumed diat small clusters of (C02)io were populated with 50% of the total 
release cloud density. The soluticm of the model is shown in Figure 28 where die density time 
histories of the important species CO 2 and (H}D) are presented. Examination of Figure 28o shows 
that a secondary peak in CO 2 density is achieved due to the breakup of large clusters into smaller 
ones. This process quickly peculates (€ 02)2 duste*^ ^ a pedc at about 100 seconds whidi then 
leads to the £ut introduction of single carbon dioxide molecules into the already relaxing diemical 
system. This creates the secondary peak in C02 density and, therefore, in 0(^Z>) as well. Figure 
2 %b clearly shows the formation of a hump in the oxygen decay curve. In addition, die maximum 
density of oxygen is reduced frcrni the previous models' results and allows an estimate of peak 
6300A enhancement to be 78 Rayleighs, only rfoout S percent above the observed maximum. On 
the basis of these results, this cluster model is a feasible description of the chemical processes that 
led to the unexpected hump and peak emissi<m observed in the Red Air I hi^ altitude 6300A 
enhancement profile. 

9.4. Sublimation of Frozen CO 3 

The Red Air 1 release conditions could have fi^ozen a large amount of the CO 2 during die 
release of the gas instead of clustering it. This suspicion arises from the feet that during the rocket 
flight the CO 2 canister was held at a temperature just above the critical temperature. Joule- 
Thonqison cooling and directed acceleration effects fh>m the fest discharge of hi^ pressure gas 
leads one to speculate that as die gas was leaving the canister it froze, seeding the release area with 
small chunks of frozen CO 2 The subsequent sublimation of the ice could prolong the introduction 
of €02 into the system and thus 0(^Z>) productiem. This effect was modeled by assuming that die 
frozen CO 2 particles were small and would sublime at a constant rate over 300 sectmds. Freezing 
and sublimation of 50% of the CO 2 is modeled by including a C02(fiozen) species and a 
sublimation reaction: 

C02(fro2en)->€02 (9.10) 

Figures 29 a and b show the results of this effort. The CO 2 density time history shows a delay 
in the decay of that species while the sublimation process continues. The resulting 0(^Z)) profile 
exhibits a clearly defined hunqi shortly after the sublimation is complete. The peak oihancement 
over the background frmn this model is estimated to be 150 Raylei^, about a frmtor of two too 
high This can be attributed to not knowing the proper way to model the sublimaticm of frozen 
release particles. The peak oihancement of the emissions is directly related to the peak of CO 2 
density. These results show freezing-sublimation of CO 2 can explain the observed 

emissions. 
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MODEL mb. CLUSTER CHEMISTRY 



TIME (SEC) 

Figure 28. (a) CO 2 cluster density time histoiy for high altitude release, small cluster removal 
model: (CX^^ + M (C02)n.2 + (C02)2 + M. 50% of total CO 2 released as (C02)2o- (b) 

0(^Z7) density time histoiy for high altitude release, small cluster removal model: (CX) 2 )ii + M -> 
( 002 ) 0-2 ((^2 released as (C 02 )io- 
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MODEL IV, SUBLIMATION CHEMISTRY 



Figure 29. (a) COj density time histoiy for high altitude release, sublimation model, (b) 0(*D) 
density time histoiy for high altitude release, sublimation model. 
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10. DISCUSSION OF COs RELEASES 

COj-iono^heric diemistry has been examined through comparison of field measurement and 
simulation results of die Red Air ICO 3 releases. Simulations that use standard COrionospheric 
reactions predict more 6300A airglow dian was observed. The simulated aiiglow is too large by a 
&ctor of fiMir fiir the upper release and a fiKtor of five for die lower release. The upper release 
may have suffered more removal of CO 3 by fireezing or clustering due to the less dc^ ambient 
gases. Hie discrepancy may be due to error in the published charge exdiange rate of O*^ widi 
CO]. The conditions were approximately diermal in the experiment with CO] heating to lOOOK. 
Due to die lack of laboratory measurements at higher tonperatures, the CO] vibrational-state 
diemistry used in this study may not represent CO]-ionospheric chemistry well. The discrqiancy 
between simulation and oqieriment also could be due to misrqiresentatitMi of the actual ambient 
conditions. Increase of the electrcm density at die upper release affects the 6300A airglow 
enhancement prc^rtionally. 

Not only do the observed and predicted airglow magnitudes differ, but also the studies of the 
airglow intensity versus time curves differ. The 6300A airglow observed at the positicm of the 
lower release is well modeled by addirtg the modeled aiiglow from upper and lower release. Two 
rodcets fiir the Red Air n releases were used to ^arate the clouds and avoid the airglow oveii^. 

The intensity-time curve of the 6300A emission frmn the upper release also has a complex 
shape. Modeling indicates that die lower release did not expand up to the upper release. The 
conqilex intensity-time curve i^ipears to be caused by other physi^ processes. Chemical and 
diermodynamic medianisms were examined as eiqilanations fiir the bimodal airglow intoisity 
curve. Explicit vibration chemistry was included but it was determined that NLTE effects could 
not produce the observed "hun^i" in the 6300A signature. It was shown that perh^s the fireezing- 
subltmation or clustering-evaporaticHi of CO] could produce the second hump of the airglow curve 
in the upper release. A simulation widi CO] clustering and evaporation produced the best 6300A 
intensity-time i ve sluqie. Significant clustering of the CO] probably occurred in the Red Air I 
release. 

The fieezing or clustering of CO] in future releases may be alleviated by heating the dispenser 
hotter than die CO] critical temperature of304K prior to release to overcome the si gnificant Joule- 
Thonqison cooling. Also, baffling the release chamber will reduce thermal cooling due to directed 
acceleration of the gas. 

Examinati<» of the ambient aiiglow signatures predicted with two sets of airglow chonistry 
(Solomon andAbreu [1988] and Unk and Cogger [1989]) suggests that atomic oxygen q uenching 
of OCD) could be inqiortant. AdditionaUy, die best match of the upper and lower release airglow 
intensity ratio was obtained using a non-zero quenching rate (at least 4x10*1^ cm^ s'^). 

Because the processes in the release clou^ included release cloud mixing and fieezing or 
clustering, definitive assessments of CO] and 0(^D) diemistry could not be made. Red Air n may 
be imne definitive in answering these questimis. 

11. HF HEATCR INDUCED CHEMISTRY IN THE UPPER ATMOSPHERE 

The chemical effects within an HF (hi^-fiequency) radar beam were studied through 
modeling of the oieigy deposition-emission processes and detailed induced chonistry. This effort 
was done in support of assessing the environmental impact of the HF Active Auroral Research 
Program Imiospheric Researdi Instrument (HAARPI^ HF facility. The HAARP radar will be 
used to heat the iimosphere to examine instabitities and other ionospheric processes. Ionospheric 
heating eiqierinients have been used primarily to study plasma dynamics and instabilities [e.g., 
Carlson and Duncan, 1971; Dubois etal., 1990]. The publish^ literature of chemical effects 







widiiii HF heater beams fi)cused primarily on F region airglow production etai, 1970; 

Stpleretai, \91A\ Bernhardt etai, 1989; and references therein]. Others have examined energy 
balance and election temperatures within HF beams [Farley, 1963; Thomson, 1970; Meltz and 
LeLevter, 1970; Meltz et al., 1974; Mantas et al., 1981; Bernhardt and Duncan, 1987; Barr and 
Stubbe, 1991; Stocker et al., 1992]. Most of the studies cited have examined F r^on enmgy 
dqpositioas. 

The set of chemical reactions used in our study of the ionospheric heating erqjcriments is 
for E and D region models. This will allow us to examine the magnitude of the HF beam 
nwtififtariftn of the natural chemical system in regums of the upper atmosphere. In particular, 
energy enhanced chemistry can produce HOx and NO^ species, whidi may subsequently be 
transpwted to lower altitudes. These species can adversely participate in the ozone chemistry of 
the stratosphere and mesosphere. 

An initial precursory effort to model the HF heating nuxlification was performed to discover 
foe possibility of a null effect on NO^ arul HO^ densities. These results will provide the 
foundation to determine the best avenue for dewloping higher-order nuxlels of the modification. 

Recent attention has been given to the effect of the energetic particle events on ozone layer 
concentration [Jackman etal., 1980; Rusch etai, 1981; Orsini and Frederick, 1982; Frederick 
and Orsini, 1982; Solomon etai, 1982; Callis etai, 1991]. The energetic particles ionize air 
molecules. The ionization and recombination chemistry produces odd nitrogen species. Important 
odd nitrogen species are N, NO, NO2, NO3, HNO3, HNO4. It was determined by Crutzen [1970] 
and Johnston [1971] that odd nitrogen plays a role in determining ozone (O 3 ) densities within the 
stratosphere and mesosphere. The amount of ozone present in the stratosphere is a result of a 
balance between processes that produce ozrme and processes that destroy ozone. 

12. ODD NITROGEN-OZONE CHEMISTRY 

The ozone destruction ability of CFCs (Chlorofluorocarbons) has been well publicized. 

Ozone is primarily created through the three-body association reaction: 

O+O 2 +M-> 03 +M (12.1) 

Atomic oxygen is created through photodissociation of O 2 by solar EUV in the upper atmosphere. 
Two important reactions in odd nitrogen-ozcme chemistry are 

NO+O 3 -►NO 2 +O 2 ( 12 . 2 ) 

NO 2 +O->N0 + 02 (12.3) 

Reactions (12.2) and (12.3) have a combined net result of: 

O 3 +O- >202 

Tlie loss of atomic oxygen reduces the rate of O3 production via reaction (12.1). The odd nitrogen 
molecutes are not destroyed in the reactions (12.2) arul (12.3). NO arul NO 2 can participate in O 3 
and O destruction repeatedly. This is the nature of a catalytic reaction. The above chemistry is a 
sinqrlification of NOx diemistry but demonstrates the main effect of increased odd nitrogen levels 
in tte stratospheric ozone layer. 

The lai]^st source of odd nitrogen in the upper atmosphere is the reaction 
N20+0(‘D)->NO+NO (12.4) 
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This source of NO produces roughly 4.SxlO^^ NO molecules per year in the stratosphere and 
mesosphere. 

Hi^-eoeigy particle events are also in^xirtant sources of odd nitrogen molecules. Enngetic 
particles generated in solar flares produce "particle events" in die upper atmosphere. The fiist 
protons in solar proton events (SPE) strike the atmosphere and produce electron*ion pairs. For 
evuy 35 eV of energy deposited in the atmosphere one electron-ion pair is produced. Other 
eaagetic particle events, such as relativistic electron precipitation (REP), and galactic cosmic rays 
(GCR), ionize air with the same efficiency, hnpact ionization reactions are primarily 


N 2 +e->N 2 +e + e (12.5) 

N 2 +e-►N'"+N + e+e (12.6) 

O 2 +C->02'^c + e (12.7) 

02 +e->0^+0 + e + e (12.8) 

0+e-►O'^+e+e (12.9) 

Ionization of air by electron or proton impact eventually produces odd nitrogen species through 
recombination of an icm-electrcm pair 

ffO*+e ->N(^D) + 0 (12.10) 

Nj-i-e-►N(^D) + N x2.1l) 

and NO production through 

n(^d) + 02-->N0 + 0 (12.12) 


Careful accounting of odd nitrogen production suggests that 1 .3 to 1.6 odd nitrogen molecules are 
produced per ionization [Rusch et ai, 1981]. The number used most often in the literature is about 
1.3. Additional sources of odd nitrogen are icNiization of air by meteors and photochemical 
production in the thermosphere The long-lived NO and N(^S) created in the thermosphere are 
subject to vertical transport. The odd nitrogai produced above 100 km altitude is transported 
downward to mesospheric (50-80 km) and possibly stratospheric (15-50 km) altitudes [Solomon et 
al., 1982]. Odd nitiogra molecules, which are created above 100 km and survive the transport to 
die lower mesosphere and stratosphere, can participate in ozone destruction. 

The only loss mechanism for NO and N during the downward transport to the stratosphere is 
die {diotodissociation of NO, coupled with a cannibalistic reaction of NO and N 

NO + hv->N+0 (12.13) 

N+NO->N2+0 (12.14) 

The photodissodation is needed to maintain N densities because N also reacts with O 2 to produce 
NO (reaction (12.12)). 
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During the long, daHc polar winters the loss rate of odd nitrogoi is very small. The 
probability that NO reaches the strat<»phere increases. Also, the vertical motion is generally 
downward in the upper atmosphere thrwgh the winter season. Solomon etal. [1982] examhied 
hi^ latitude odd nitrogen promotion and transport. They concluded that auroral activity during 
winter months can affect stratospheric odd nitrogen levds. They also conclude diat sununo’ 
hemisphere odd nitrogen densities in the stratosphere are not appreciably enhanced by auroral 
activity, because the N and NO densities ate brou^t into equilibrium Arou^ photodissodation 
and subsequent mutual destructirm during downward transport. These conclusions are suj^rted 
by rqmtts of enhanced NO densities in the winter polar r^on [Rusch and Barth, 1975]. 

A comparison of natural odd nitrogen sources with the nuclear test series odd nittx)gea 
production provides necessary context in understanding the magnitude of possible perturbatioos 
caused by hiunan activity. Jackman and others [1980] compile a list of tte armual production of 
odd nitrogen fixmi various sources. Tables 9 and 10 contain the results of the crmpilation. The 
nuclear tests produced nearly as much odd nitrogen as the largest natural source. The tests were 
probably resprmsible for a 5% reduction of global ozone layer content. 


Table 9. Stratospheric NO sources [Jackman et al ., 1980]. 


Source 

NO nwleculestyr 
solar maximum 

NO moleculestyr 
solar minimum 

N2QKX^D) 

4.5x1034 

4.5x1034 

REPS 

2.7x1031 

2.7x1030 

SPEs 

2.5x10^3 

2.5x1033 

OCRs 

2.7x1033 

3.7x1033 

Nuclear testing 1961-62 

2.2x1034 

2.2x1034 

Total without Nuclear tests 

5.0x1034 

5.0x103'» 


Table 10. 

Mesospheric and thermo^heric NO sources [Jackman et at., 1980]. 

Source 

NO moleculestyr 

NO moleculestyr 


solar maximum 

solar minimum 

N20+0(1D) 

5.3x1032 

5.3x1032 

REPs 

1.4x1034 

1.4x1033 

SPEs 

3.9x1033 

3.9x1033 

Thermospheric 

1.5x103-1 

3.7x1033 

Nfeteors 

6.3x1032 

6.2x1032 

Total 

3.4x103-1 

1.0x1034 
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13. ENERGY DEPOSITION MODEL 

There are two known mechanisms for the omversion of HF field energy to particle eneisy: 
ohmic and anomalous. In an underdense itmosphere (that is, the HF emission firequency is much 
less than the plasma fiequency) the absorption of the HF beam takes place through ohmic heating. 
Beam energy is absorbed as accelerating electrons collide with ions a^ neutral constituents. The 
directed energy of the accelerating electrons is dMrmalized in the collisicms with neutrals and ion as 
well as self-collisions. On a much slower time scale, the ions and neutrals are heated fiirough 
elastic and inelastic collisions with the electmis. Additionally, anomalous heating occurs when die 
plasma fiequency approaches the HF enhssion fiequency. Anomalous heating results fixxn time- 
varying plasma instabilities and occurs in a very narrow altitude range. 

Typical HF radar specifications are listed in Table 11. Ambient backgrounds for winter 
nighttime are used for tte investigations. 

Table 11. Typical HF radar qiecifications. 

Total Beam Power 3MW (2x10^5 eV/s) 

Beam width 10° to 3° 

Efieedve Beam Power 0.7 to 3GW 

Frequency Range 10 to 3 MHz 

Yeariy operation 320 hrs n.lSxIO^ s) 

Yearly emissioa 3.5x10“ W/year (2x10^^ eV/year) 


Ambient cemditions for the modeling of HF effects were obtained firom MSIS90 neutral 
atmo^ihere [Hedin^ 1990], the E region model of Rasmussen et al. [1989]. Electron profiles of 
Ogawa and Shimazah [1975] are used below 90 km. The geophysical parameters of average solar 
activity and moderate magnetic storm activity are used in obtainiiig the neutral atmospheric and 
ionospheric profiles. 

Ohmic heating of the electron gas is the assumed heating mechanism for energy dqxrsition in 
the D and E r^(» of the ionosphere. The literature has many examples of deposition (Ulculadons 
[Farley, 1963; Thomson, 1970, Meltz and LeLevier, 1910; Showen and BehrAe, 191%; Perkins 
andRohle, \919;Mantas etal, 1981]. The rqiproach of Perkins ondRoble [1979] is followed 
and extended for the calculations presented here. For the polar ionosphere the electron energy 
equation can be eiqiressed as 

-I - i., - ^ (.3.1) 


where t is time, z is altitude, Ug is electron density, K is Boltzmarm's constant, Tg is electrcm 
tenqierature, Ug is electron velocity, “kg is thermal conductivity, Qjjp is the HF energy dqxisition 
term, Lgi represents all elastic collisional energy losses, L;„gi represents all inelastic collisional 
energy losses, and represents all ionizing collisional energy losses. The ohmic HF heating 
term is; 


Qhf - 



(13.2) 


where F is the energy flux of the beam, c is velocity of light, v^ is the plasma frequency, v^ is the 
beam fiequency, and Vg is the total electron collision firequency (that is, combined electron-neutral 
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collision fiequency, Vg„, and electron-ion collision frequency, v^,). The energy flux, F, is 
with altitude as beam energy is absorbed. A kilometer below the altitude M/bere 
approaches Vp significant energy dqiosition occurs dim to anomalous heating. The anomalous 
heating can produce non-fiiennal electrons with energies of 2 to 10 eV [see, Bernhardt et al. 
[1989] and references dieiein]. The anomalous beating within the HF beam is observed in the F 
region. A large percentage of die beam energy is absmted in the D and E region wiiere ohmic 
heating dominates. 

The heated electron gas is cooled through the transfer of energy to the heavy particles by 
elastic and inelastic collisions. The elastic coUisional heat loss is writtoi 


iw = z«. ^ v„3Ji:(r. - j;)+ ^ v„3i(:(r. - r,) 

i ^ n ^ 


m. 


(13.3) 


bielasdc heating is considerably more complex. Perkins and Roble [1979] include excitation of 
rotational and vibratkmal modes of O2 and N2, fine structure levels of O, and the metastable 
species 0(^D). The vibrational and rotadmial eiwrgy of O 2 and N 2 is quickly transferred to the 
^cient radiating mdecules CO2 and NO. Electronically excited molecules radiate or quench 
collisionally. 

Perkins and Roble [1979] noted that the coUisional cooling of the electr<m gas is not sufficient 
to keqi die electron temperature fiom "running away" from the neutral temperature. The "run¬ 
away" electron temperature is capped at higher eneigies through inelastic coUisions that have larger 
energy loss per collision (i.e., impact excitation of higher electronic states and impact ionization). 
Additional excited states are included in the inelastic loss term. An impact-ionization cooling term 
isalsoadded. The vibrational and rotatirmallossrates of5<rAunikora/Migy[1978] are used. The 
excited states included in die model are listed in Table 12. The electron impact excitation cross 
sections were obtained firom Banks andKockarts [1973]. The energy loss terms for the electron 
gas duough inelastic collisions can be written in the fomi 

^»€l=Einel^weM (^ 31 ) 

vdiere £/„«/ is the energy of the excitation of the mutral species, ki„gi is the rate of the excitation of 
die neut^ species, and n„ is the density of the neutral species. 

The dectron density is can be sensitive to energy dqiosition rate. As the electron density 
dianges so does die energy dqiosition rate. A self-ccmsistent determinatimi of the electrrxi density 
and the energy dqiosition is necessary for a ctmqilete modeUng effort. For example, within a low 
power HF beam, electron energy may not be sufficient to ionize but the elevated electron 
tenqierature will decrease the ion-electron rectHubination rate. This increases the E region plasma 
density. As the HF beam power increases and the electron energy with it, impact ionization may 
beemne important. The electron impact itxiization cr(»s sectiims were obtained from Banks and 
Kockarts [1973]. The loss terms for impact icmization have the form 

Li^^Eipa„ngn„ (13.4) 

idiere Ejp is the iemization potential and ag„ is the impact ionization rate. 
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Table 12. Electronically excited species contained in the deposition model. 


Gas 

State 

Energy (eV) 

Lifetime (s) 

N2 

a3i„+ 

6.14 

1.9 


B^n- 

c’nS 

a^n- 

b*n5 

7.30 

8.9x10-6 


11.03 

3.6x10-6 


9.10 

1.5x10-^ 


12.85 

10-"^ 

02 

a^Af 

.98 

3880 


1.64 

11.8 



4.5 

10-"^ 



8.3 

lO'"^ 

0 

Id 

l.% 

148 



4.17 

.8 



9.53 

1.8x109 


5s 

9.15 

6x10-4 


14. RESULTS OF HF HEATING MODEL 

The two electron density profiles are used in modeling the HF heating effects (Figure 30). 

Both profiles are fiom polar winter. For simplicity the initial electron temperature was set equal to 
die noitral temperature at all altitudes. The energy depositimi model described in previous sections 
calculates magnitudes of energy deposititm, electron temperature and cooling terms for HF heater 
beams. 

Energy deposition, beam flux, electron tenqieratuie, and ionizadmi rates for the two electron 
density profiles are shown in Figure 31. The energy dqiosition drc^ off sharply below 75 km for 
both cases due to the absence of electrons below this altitude. The electrons attach to molecules 
readily dirougfa diree-body attachment. During daylight, photodetachment increases the electron 
density between 60 and 75 km altitude. However, bram attenuation below 75 km during daylight 
will be mwiimiTied as electron attachm^t is accelerated by die elevated electrcm temperature in the 
HF beam. When the beam is on, the free electrcm density within the beam will be reduced below 
75 km attitude. Three-body attachment is not important above 85 km. 

The denser ionosphere within die auroral o\^ reflects die 2.8 MHz beam so the flux drops to 
zero (Figure 31c). Greater energy deposition rates occur in die icmosphere of the auroral r^on. 
The beam attenuation is obvious in all examples except the 8MHz beam in the less dense 
ionosphne (Figure 31c,r0- 

Electron tenqieratures of the D and E regicxis calculated fiom the balance of cooling and 
heatiiig terms are larger than the temperature enhancements calculated in the Perkins and Roble 
[1979] model of Arecibo HF heating experimoits. The same total onissicm energy was used 
(3MW) but the beam width is narrower that the previous study providing for larger energy fluxes. 
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Figttfe30. Electron density profiles in and out of the auroral oval. The conditions are for 
nighttime winter with moderate auroral activity. 


A very hi^ electron temperature results in ionization and dissociation of the neutral 
atmo^hete diroug)i electrcm impact. The temperatures of .8eV at 90 km within HF beam are still 
not hi^ enou^ for ionization or dissociation to become important. Even though the local energy 
deposition of the HF beam is comparable to an intense aurora, the heated electrons do ionize 
efficiently (Figure 32g,/t). Most of the energy goes into vibrational and rotational excitation of O 2 
and N 2 . Hie ratio of the rates of ionization to energy dqposition determined from the energy 
dqxrsition model predicts one ion-electrmi pair produced per 10^^ eV energy deposited. The 
energy dqiosititm from high-energy particle events (>S00 eV), such as in auroras, SPEs, GCEs, 
RPEs, prices ion-electron pairs efficioitly. For such events one ion-electron pair is created per 
35 eV deposited. 

The in^MCt of HF heating on the odd nitrogen chemistry can be estimated from the ion- 
electron pairs produced within the beam. Using the number of 1.3 odd nitrogen molecules per ion- 
electron pair, the total beam emission rate and 10^^ eV per ion-electron pair, 2x10^3 molecules/s 
are produced in the beam. We can now cmnpare these rates to the natur^ production mechanisms 
of odd nitrogen listed in Tables 9 and 10. For tlm typical operation of an HF radar (Table 11) an 
armual operation time of 10^ s is expected. Therefore, 2x10^^ odd nitrogen molecules could be 
produced during a year of HF operation. This is 14 orders of magnitude less than natural 
processes. 

The possibility anomalous heating within the HF beam was mention previously. Anomalous 
heating can produce dectrons with 2 to 10 eV orergy within the F r^on of the ionosphere. This F 
r^on heating is only a limited percentage of the beam energy. Much of the beam energy will be 
dqxrsited in foe D and E region vfoere dimic heating dominates. The higher energy electrons of 
anomalous heating are still less efficient in ionizing air than particle precipiation events. The ratio 
of one ion-electron pair per 35 eV is used for particle energies above 500 eV [Banks and Kocharts, 
1973]. 


57 





E 

jt 

I 

I 


200 

ISO 

160 

140 

120 

100 

SO 

60 

40 

20 

s 


N. Insid* Ovoi . 

2.auH2 — 

200 

1B0 

160 

_ 140 

J. '20 

€ 100 
3 

Outside Oval . 


3 ftO 

< 

■-— 2 SMHz 


60 

* 


40 

■ 

'.—.-__ __ 

20 

0 

• 


0» 10J 10* 10* 10* 101 10* 10* 


10 » 101 10 * 10 » 10 * 101 10 * 10 » 


E 

JK 

V 


E 

JK 

I 



T,(K) 



200 
180 
160F 
140 
1 20 
100 
80 
60 
40 

20 F 
0 


Oulside Ovoi 


(H) 




'8MH^ 


2 8MHz 


10-»1 0-41 0-110-* 10-110-*'0-110-110-' IC'^ 
lonizotion Role (cm-’s"') 


Figure 31. Results ofdepositioa model are for beam frequencies of 2.8 and 8 MHz. Local energy 
dqiosition for (a) inside the auroral region and (b) outside the auroral region. Energy flux of HF 
beiuns for (c) inside the auroral region and (d) outside the auroral region. Electron temperatures 
fc' (e) inside auroral region and (f) outside auroral region. Ionization profiles for (g) inside the 
auiora zone and (h) outside the aurora zone. 
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15. DETAILED CHEMISTRY OF A PARCEL OF AIR IN AN HF BEAM 


The tiine<dependent detailed chemistry of a parcel of air within the HF beam can be studied 
using infoimation obtained from the eneigy deposition model described in the previous sections. 
The detailed chemistry model is used to examine the efficiency of odd nitrogen molecule creation 
within an HF heater beam. The energy deposited within the HF beam will not produce ion-electron 
pairs. However, c^r reactions could be driven by the elevated electron temperatures Are there 
other avenues of odd nitrogen production within the beam? This is the question to be answered in 
this section. 

The calculated results of the energy deposition model of the HF beam are used to drive the 
diemistryofaLagrangian parcel ofair moving through the beam. At 90 km the HF beam, where 
the dqrosition peaks, is S km in diameter. The wind velocities vary between 0 and 60 m/s with 35 
m/s average at 90 km altitude. The electron temperature and energy dqx>sition histones for the 
Lagrangian parcel are shown in Figure 32. The initial densities for the detailed chemistry study are 
obtained from MSIS90 [Hedin, 1991] and a model minor species mixing ratios [Anderson et ai, 
1986]. Ionospheric constituents and ionization source terms are from the Rasmussen et al. [1989]. 
A fully time-dqrendent chemistry model and an extensive reaction database is used to solve a 
system of ordinary differential equations for 152 atmospheric species and 1844 reactions. The 
chemistry is defined by a comprdiensive set of reactions [Kennealy et al, 1989; Laridn, personal 
communication] for the regions of the HF radar energy deposition. The conditions in the air parcel 
are allowed to arrive at steady state before the energy dqiosition is increased. Two cases are 
examined. The first case represents a high-energy particle precipitation event that varies in time 
and the second case represents energy deposition in a moving parcel of air within an HF heater 
beam. The energy deposition profile is the same for both cases (Figure 32a). Only the ionization 
efficiencies are different. 



Time (s) Time (s) 


Figure 32. (a) Energy deposition rates and (b) electron temperature seen by a parcel of air moving 
through the HF heater beam at 30 m/s. These values are us^ to drive a detailed chemistry morfel, 
first, within an HF heater beam and, second, during a high-energy particle precipitation event. 
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For die first case of the high-energy particle precipitation, the energy dqrositioa curve of 
11. la is used to drive ionization reactions wdrose rates are defined by the "magic number", 35 eV 
per one ion-electron pair. The resulting chonical changes during die high-energy particle 
precipitation event are plrmed in Figure 33 with dotted lines. The increase in plaana density is 
large during the precipitation event but the plasma relaxes to previous densities with some 
differences in the pre- and post-event conqxisition (that is, less O 2 ''' after the event). The odd 
nitrogen species, NO, N, and NO 2 , all increase dramatically. Only N decreases after the event. 

NO and NO 2 remain at the elevated densities due to the absrace of destructimi processes during 
polar winter. 

The second case examined was a parcel of air moving through the HF beam at 30 m/s. The 
electron temperature curve in Figure 32a was used to drive the chemistry. The ele^ron 
temperature history (Figure 32b) was derived from the energy dqiosition curve show in Figure 
32a. The neutral and ion temperature histories were also included in the model. Botii tonperatures 
only increase about 1° K during the parcel transit of the beam. 

Within the HF beam, the elevated electron temperature slows tiie dissociative recombination 
rates. This increases the plasma density (Figure 33a). The increase in plasma density is not due to 
increased ionization levels. After leaving the beam, the plasma within the parcel returns to ambient 
dmsity aixl composition. The odd nitrogen levels are imchanged throu^out the transit of the beam 
(Figure 33b). There are some composition changes vthen the parcel is in the beam, such as 
increased metastable species, but these return to ambient levels upon beam exit or turn off. There 
is no observable odd nitrogoi production. An NO production between tiie two cases of energy 
deposition is more likely due to non-ftiermal characteristics of electrons in electric fields [Carron, 
1992]. 

The small region within the HF heater beam has local energy depositimi rates that are 
comparable to intoise auroras. This permits well determined modification experiments to be 
performed on the ionosphere. In particular, energy deposition and excited species chemistry might 
be examined through carefully plaiuied optical experiments. The studies may illucidate the 
chemical time constants within the E and F regions. 
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Figure 33. The chemical changes are shown for two different cases. The energy deposition is the 
same in both cases but one case represents HF heating effects (solid lines) and the other represents 
high-energy particle precipitation effects (dotted lines), (a) The main plasma constiments and (b) 
the odd nitrogen molecules within the air parcel are shown for the two cases. 
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16. DISCUSSION OF HF HEATING 


We have conqjaied atmospheric oiergy dqx^ition of nature to that of an HF radar to 
determine the magnitude of loQg term chemical effects of HF operations in the atmosphere. The 
total energy emitted by an HF radar is many orders of magnitude smaller than the mesopheric and 
thermospheric energy deposition by natural processes. ^ gy deposition, whi^ leads to 

ionization and odd nitrogen production, is also orders oi larger that the yearly energy 

emission of an HF ftcility. By assumirig the most efiBcient production of odd nitrogen per energy 
deposited we show that ^ HF radar has essentially no efife^ on the global chemical halanc<^ 
However, the HF radar can be a perturber of local quantities for r^ons of lO’s of kilometers. The 
perturbirig ability may be an important tool for investigating D and E region chemistry. An 
important next step frcmi this study is to include non-maxwellian electron distributions [Carron, 
1992] to estimate foe NOx production efficiency within HF beams over a range of conditions and 
modes of (^ration. 
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